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The Debye theory of polar molecules leads to the conclusion that the anomalous 
dispersion due to dipoles which has been observed at high frequencies should also be 
found at low frequencies. Supercooled glycols are viscous enough at low temperatures 
to show dispersion at low frequencies. Real and imaginary dielectric constant and 
power factor are calculated and their inter-relationships are discussed. Better than 
qualitative agreement with the Debye theory is observed in the pure pentandiol, but 
not in the solution measurement. 


T IS well known that in applying the Debye equations for determining the 

electric moment, the results obtained for pure liquids, especially those 
ordinarily considered to be associated, do not agree with those found for 
dilute solutions or gases. Although the simplifying assumptions of the Debye 
theory are no longer to be regarded as valid for associated liquids or for solids, 
such a treatment appears qualitatively correct and is useful as a means of ex- 
plaining the effect of temperature and frequency upon dielectric properties. 
For the purpose of relating dielectric properties to chemical and physical 
structure, measurements on pure liquids are quite as important as dilute so- 
lution measurements. 

A study of the anomalous dispersion, i.e., change of dielectric constant 
and dielectric absorption with frequency, provides a useful test of the Debye 
theory. This change of dielectric constant with frequency in the electric 
spectrum is usually called anomalous dispersion to distinguish it from the 
normal dispersion observed at optical frequencies. Anomalous dispersion is 
observed for highly damped polarizations, as for example for dipole polariza- 
tion, where the time of relaxation is determined primarily by viscous forces. 
Normal dispersion is found for lightly damped polarizations such as elec- 
tronic polarization where the inertia of the charges is the determining factor. 

Most of the measurements of the variation of dielectric constant with fre- 
quency which have been‘made on pure materials have been made at high fre- 
quencies. Mizushima’s! measurements on a series of alcohols, for example, 
were at 0.57, 3, 9, and 50 meters. Johnstone and Williams? measured solu- 


1S. Mizushima, Bull. Chem. Soc. Japan 1, 47, 83, 115, 143, 163 (1926). 
2 J. H. L. Johnstone and J. W. Williams, Phys. Rev. 34, 1483 (1929). 
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tions of nitrobenzene and para-dichlorobenzene in a viscous oil at several fre- 
quencies between 34 and 154 meters. Goldammer and Sack,’ working at 52.8, 
77.5 and 139 cm, showed a change of dielectric constant with frequency for 
solutions of some alcohols ina mineral oil. Of these investigators only Mizu- 
shima has given conductance data from which the absorption accompanying 
the anomalous dispersion can be determined. There have been numerous 
measurements‘ on rosin, insulating oils and mixtures of the two, but the 
structure of these materials has usually not been known. 




















TABLE I. 
Substance mr. Np?” ad 
Propylene glycol 65°-66° at 9 mm 1.4323 
Trimethylene glycol 87°-88° at 9 mm 1.4396 
2-Methyl pentandiol 2.4 92°-93° at 7 mm 1.4316 








There is no reason from the theory for believing that the same anomalous 
dispersion shown for example by Mizushima for the alcohols, would not also 
be found at low frequencies. Most materials do not show this dispersion at 
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Fig. 1. Dielectric constant and power factor vs. temperature for frequencies from 
1 to 100 kilocycles for propylene glycol. 





low frequencies due to the fact that their viscosity is not sufficiently high. 
However, such polar liquids as may be supercooled to form glasses usually 


3 R. Goldammer and H. Sack, Phys. Zeits. 31, 224 (1930). 

4 D. W. Kitchen and H. Mueller, Phys. Rev. 32, 979 (1928); E. Kirch and W. Riebel, Arch. 
f. Electrot. 24, 355, 553 (1930); L. S. Ornstein and G. Willemense, Zeits. f. tech. Physik 11, 
345 (1930); H. H. Race, Phys. Rev. 37, 430 (1931). 
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do show anomalous dispersion. The glycols provide a series of substances 
which become sufficiently viscous at low temperatures to show this dispersion 
at low frequencies and this paper describes the results of measurements on 
three glycols. The characteristics of the glycols used are given in Table I. 
These materials were dried and then distilled at low pressure. 

The measurements were made with a capacitance bridge described by 
Shackelton and Ferguson® and were direct capacitance and conductance 
measurements. They were usually made with a descending order of tempera- 
tures but sufficient time was allowed after changing the temperature for the 
dielectric constant and conductance to become constant. That this was so is 
indicated by the fact that points taken with ascending order of temperatures 
fall on the same curves as those taken with descending temperatures. 
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Fig. 2. Dielectric constant and power factor vs. temperature for frequencies from 
1 to 100 kilocycles for trimethylene glycol. 


In Figs. 1 to 4 are shown curves of dielectric constant and power factor 
for the three glycols and the solution of pentandiol in the non-polar solvent di- 
oxane. The points are the observed values. It was not possible to observe any 
dispersion for very dilute solutions of these glycols although, as can be seen 
from Fig. 4, essentially the same behavior is observed for a concentrated solu- 
tion as for the pure glycol, it merely being displaced to a lower temperature as 
required by the theory. Bridge measurements can not be made at frequencies 
such as would be necessary to observe dispersion for the dilute solutions and 
high frequency conductance measurements by other methods are more diffi- 
cult to make and less dependable. 


' W. J. Shackelton and J. G. Ferguson, Bell System Tech. J. 7, 70 (1928). 
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It is apparent from the curves that the different glycols, and even the 
solution, all show the same general behavior. The curves are, however, dis- 
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Fig. 3. Dielectric constant and power factor vs. temperature for frequencies from 1 to 100 
kilocycles for 2-methyl pentandiol 2.4. 


placed along the temperature axis, the dispersion of dielectric constant first 
appearing for pentandiol, Fig. 3, at —35° while for trimethylene glycol, Fig. 
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Fig. 4. Dielectric constant and power factor vs. temperature for frequencies from 1 to 100 
kilocycles for 0.83 MF solution of 2-methyl pentandiol 2.4 in dioxane. 


2, it starts at —60°. At the higher temperatures the dielectric constant of 
each of these materials is independent of frequency in the range below 100 
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kilocycles and is increasing rapidly with decreasing temperature. As the tem- 
perature is lowered the viscosity increases and the time of relaxation of the 
dipole, which is given by the equation 


dra*n 
T= ; (1) 
kT 

where 7 is the viscosity, a the molecular radius and 7° the absolute tempera- 
ture, increases. When this becomes of the same order as the period of the ap- 
plied field a decrease of dielectric constant is observed. At this temperature 
also the power factor starts to increase and it reaches a maximum near the 
temperature where the change of dielectric constant with frequency is great- 
est. 

The general equation for the variation of dielectric constant with fre- 
quency as given by Debye is 


€0 : € 
- + wr 
(€9 + 2) * 2) 
€ = 
. 1 
+ wr 
(€o + 2) (€,, + 2) 


where € is the static dielectric constant or the value at a low enough fre- 
quency that dipole orientation is constant, and ¢, is the high-frequency di- 
electric constant where the effect of dipole orientation has disappeared. The 
equation for the dielectric constant may also be written € =e’ —ie’’. The equa- 
tions for e’ and e’’ in terms of €, €,, and 7 are 


€y —— c. 
= + 
€0 +- 2 ; 
1+ wr? 
€é. + 2 
€0 4 2 
{wT 
é.7T 2 
e’’ = (€9 — €,) 


€y + Z 2 
1+ | | wr? 
€é. + 2 


For the case of a condenser shunted by a resistance, e’ = C/C» and €’’ = 
G/wCo, C being the total capacity, Cy the capacity of the empty condenser and 
G the equivalent parallel conductance. This conductance is the non-d.c. con- 
ductance. The term e’’ which is a measure of the dielectric absorption should 
be a maximum when 


€) + 2 


The curves in Figs. 1 to 4 show the maxima for each frequency and also show 
that this is displaced toward lower temperatures (greater viscosity) for low 
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frequencies. From the position of these maxima and the values of €,, and € , 7 
may be determined and thus the general curve of variation of dielectric con- 
stant with frequency calculated. This has been done for all of the pure glycols 
and for the solution, and very good agreement between calculated and ob- 
served values is found. Table II gives the values of €,, €)9 and 7 used in cal- 
culating the theoretical curves. 


oo?’ 














TABLE IT. 
Substance Temp. €0 _ ee 
Propylene glycol —75° 59.5 3.5 2.exae * 
Trimethylene glycol —90° 66.4 : 3.0X10-5 
Pentandiol —60.4 39 2.7 1.010- 
0.83 MF solution of pentandiol in dioxane —65 35.1 2.6 2.41076 








In Fig. 5 the solid curves of e’ and e’’ are calculated from the theory for 
the pentandiol at —60.4°C and the points are the observed values. The agree- 
ment is very good. Fig. 6 shows similar caluclated curves and experimental 
points for the 0.83 MF solution of pentandiol in dioxane, and the agreement 
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Fig. 5. Curves of e’ and e’’ vs. frequency for 2-methyl pentandiol 2.4 at -60.4°C. 


is only qualitative. For the propylene and trimethylene glycols the agreement 
is also only qualitative but it is good enough to regard the Debye theory as 
being very useful in explaining the dielectric behavior of polar liquids. 

The case of the solution of pentandiol in dioxane, a nonpolar solvent, is 
particularly worth comment since it is to be expected that the best agree- 
ment between theory and experiment would be found for dilute solutions, 
where the assumptions of the Debye theory are most nearly applicable. The 
agreement here is the poorest of any. This solution, however, is not in any 
sense a dilute solution and the peculiarities of the curves of polarization vs. 
concentration for alcohols as shown by Smyth and Stoops* and for glycols as 
will be shown in a later paper show that complex associations or compound 
formations are taking place in such solutions so that the observed time of re- 
laxation is probably not that of a single molecule but an average value of 


6 C. P. Smyth and W. N. Stoops, J.A.C.S. 51, 3312 (1929). 
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that for several combinations of molecules each having different times of re- 
laxation. 

In Table III are given the values of the electric moments as determined 
from dilute solution measurements in dioxane and in one case in n-heptane. 
The third and fourth columns respectively give the values of the maxima of 
power factor and e’’ for the pure glycols and the solution of pentandiol in 
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Fig. 6. Curves of e’ and e’’ vs. frequency for 0.83 MF solution of 2- 
methyl pentandiol 2.4 in dioxane at —65°C. 


dioxane and the fifth the dielectric constant at —30°. There is no quantita- 
tive agreement between the values of dielectric constant or power factor of 
the pure glycols and the electric moment. As was pointed out earlier, it has 
been found, quite generally, to be impossible to deduce the properties of as- 
sociated liquids from dilute solution studies. 











TABLE III. 

—_ Power factor ———— 7 ; 

Substance uw X 1018 max. “= €_39° 
Propylene glycol Be 0.71 24 43 
Trimethylene glycol auae 0.77 28 46 
2-Methyl pentandiol 2.4 (dioxane) 2.9 0.71 17 31.7 
2-Methyl pentandiol (n-heptane) 2.0 
0.83 M.F. solution pentandio 0.64 12 28 





The values of the power factor maxima for the pure glycols are not very 
different from each other and may represent the range which is character- 
istic of glycols. The power factor as calculated here is G/(G?-+a*C?)!/? where G 
is the equivalent parallel conductance. In all cases the direct current con- 
ductance had decreased at temperatures where absorption was observed until 








320 A. H. WHITE AND S. O. MORGAN 


it amounted to only a few percent of the total for frequencies above 10' 
cycles. The power factor is the ratio of the in-phase current to the total cur- 
rent and is not directly significant in connection with the Debye theory. 
What is more significant is the value of e’’, the imaginary part of the complex 
dielectric constant. It is apparent from Table III that the differences in the 
value of the total dielectric constant (fifth column) are more nearly reflected 
in the value of e’’ than of power factor. 

As has been pointed out before there is some question as to the applica- 
bility of the Debye equations to such highly associated substances as the 
alcohols and glycols. Although the limitations of the theory may preclude its 
quantitative application and its use for the calculation of molecular con- 
stants from data on pure liquids, still it seems that qualitatively at least the 
rotation of dipoles may be offered as a valid explanation of the anomalous 
dispersion and absorption in substances known to contain polar molecules. 

There is the possibility that this anomalous dispersion and absorption 
which has been observed for alcohols, glycols and other polar substances may 
be due to an ionic polarization rather than to dipoles. A brief discussion of 
this possibility is given in a paper by the authors’ reporting similar measure- 
ments for rosin oils. The ionic polarization is an especially attractive explana- 
tion only at low frequencies. It can hardly apply at high radio frequencies. 
Our results, as well as those of others show however that the value of the 
maximum of absorption is essentially the same at low frequencies as at high 
although our range of frequencies is not high enough to exclude definitely the 
possibility of ionic polarization. This is the behavior predicted by the Debye 
theory but is not what would be expected if ionic polarization were respon- 
sible for the anomalous dispersion at low frequencies and dipole polarization 
at high. The constant value of the maxima of the absorption from 10* to 10° 
cycles in the case of the glycols is strong evidence that the cause is the same 
at all frequencies namely the rotation of dipoles. 

One further test of the applicability of the Debye theory to these results is 
found in the comparison of the time of relaxation as obtained from the dielec- 
tric data with the value calculated from molecular radius and viscosity by 
Eq. (1). Alternatively the molecular radius may be calculated from the di- 


TABLE IV. Viscosity and calculated molecular radius for 2-methyl pentandiol 2.4. 











Temp. n (Poises) Tr a (cm) 

—45.4 1.77 X10 6.6107 3.410-% 
—49.8 5.75 X 10° 1.810 3.210-& 
—54.9 2.5 X10 5.8107 3.01078 








electric data. Such calculations are to be considered as interesting indications 
rather than proof of the correctness of the theory. A discussion of the limita- 
tions involved in the use of the ordinary coefficient of viscosity for such calcu- 
lations is given by Williams and Oncley.* In Table IV are given the values of 


7S. O. Morgan and A. H. White, J. Franklin Inst. 
’ J. W. Williams and J. L. Oncley, Jour. of Rheology 2, 271 (1931). 





















DIELECTRIC PROPERTIES OF SOME GLYCOLS 321 





viscosity and calculated molecular radius for the pentandiol at three low 
temperatures in the range where dispersion was observed. 

Such calculations are not considered by the authors to be more than an 
approximation but it appears significant that the results are of the correct 
order of magnitude. Mizushima’s calculation of the molecular radius for the 
alcohols also agreed well with the values obtained by other means but his 
value for glycerine was too low. With the exception of glycerine the calcula- 
tions of molecular radius from dielectric dispersion data for pure substances 
have given results of the right order, which should be considered more than 
fortuitous. 

The time of relaxation of the dipole has been related to the viscosity and 
molecular radius by Debye as given in Eq. (1). This equation should not be 
expected to be rigidly applicable for associated substances such as the glycols. 
Williams and Oncley* have pointed out the difficulties involved in the use of 
the ordinary coefficient of viscosity for calculating molecular radius from dis- 
persion data. We have measured the viscosity of pentandiol in the range of 
temperatures where dispersion was observed and from these data and values 
of 7, determined from the positions of the absorption maxima, have calcu- 
lated its molecular radius. The result is a value of 0.75 X10~-° cm. This is too 
low but the result is not unexpected. It is interesting to note that the value 
falls between that of 0.35 X 107‘, calculated by Mizushima from similar data 
for glycerine, and 2.0X10-*, his value for n-propyl alcohol. 
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Development of Polar Characteristics in Insulating Oils* 


By W. N. Stoops 
Research Laboratories, Westinghouse Electric and Manufacturing Company 


(Received March 7, 1932) 


It is shown that the decomposition in transformer oils may be followed rather 
closely by measuring the dielectric constant of the oils as they deteriorate in service, 
the method being sufficiently accurate to distinguish between different oils. These data 
are compared with those obtained on oils decomposed very much more rapidly in 
laboratory tests. The change in polar characteristics is definitely correlated with 
change in chemical characteristics. That the change in dielectric constant of the oils 
is due to the formation of polar molecules is shown by a study of the variation of the 
dielectric constant and power factor of the decomposed oils with temperature and 
frequency. 


INTRODUCTION 


HE oil used in oil-filled transformers and circuit breakers is a highly 

refined product containing little or no unsaturated material. Saturated 
hydrocarbons have been shown to be nonpolar,' so that insulating oils are 
practically non-polar, that is, the molecules are electrically symmetrical. 
However, in use such oil gradually decomposes, largely due to oxidation, and 
forms various oxidation products, chiefly acids, esters, and some solid ma- 
terial or sludge, eventually becoming unfit for use. These decomposition prod- 
ucts are all polar. Hence it appeared likely that a study of the change in pro- 
perties of the oil such as dielectric constant and power factor that depend on 
polarity might be useful in following the course of the decomposition as it- 
occurs in a transformer. The following results were obtained in a study of this 
matter. 


EXPERIMENTAL METHOD 


The dielectric constant measurements were made with a capacity bridge 
previously described.? A change in dielectric constant of one part in five hun- 
dred was easily detectable. The measurements were all made at 300,000 
cycles. 

A resonance circuit recommended by the American Society for Testing 
Materials* was used to make the power factor measurements. It was modified 
by introducing an intermediate tuning circuit between the oscillator and the 
measuring circuit to remove any harmonics in the wave form. The same cell 
and method of temperature control? were used for both dielectric constant 
and power factor measurements. The power factor tests were made at two 
frequencies, 500,000 and 1,000,000 cycles. 


* Scientific paper No. 581. 

1 Smyth and Stoops, J.A.C.S. 50, 1883 (1928). 

2 Stoops, J. Phys. Chem. 35, 1704 (1931). 

3 Proc. A.S.T.M, Vol. 31 I. 877, and 1931 T.S., 642. 
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The three oils studied were first oxidized by heating a small sample of each 
in an open bottle in an oven kept at 160°C. The dielectric constant of each 
oil was measured every day or two for the duration of the test. For this pur- 
pose a 20 cc portion was withdrawn, its dielectric constant measured at 30.0° 
C and then returned to the remainder of the oil in the oven. To determine the 
importance of the oxygen in the air in the development of polar characteristics 
in the oils, glass buibs were secured in which an atmosphere of nitrogen could 
be maintained above the oil, and from which samples could be withdrawn for 
testing without introducing any air. The oils were as before all heated at 
160°C. Metallic copper is known to be an effective catalyst in the decomposi- 
tion of transformer oil. Accordingly two samples were heated under the same 
conditions obtaining in the first set of experiments, except that each bottle 
contained about two feet of bare copper wire. One of these samples was kept 
closed and the other left open during the heating. 











2 4 6 8 10 12 4 6 18 20 22 24 26 26 32 4H 3H 8B 
DAYS 


Fig. 1. Oils heated in air. 


MATERIALS 


The oils used in these tests were of different compositions. Oil No. 1 was 
prepared from a mid-continent crude and consisted largely of naphthenic 
hydrocarbons with about 20 percent of aromatic hydrocarbons and a very 
small amount of resinous material. No. 2 was a more highly refined oil from 
the same crude, being almost 100 percent naphthenic hydrocarbons. Oil No. 3 
was prepared from a Pennsylvania paraffin base crude, and was composed al- 
most entirely of paraffin hydrocarbons with no unsaturated, aromatic or 
resinous compounds. 

These three oils were chosen for study because of their different composi- 
tions and because they had been used in transformer life tests. That is, their 
decomposition in service had been studied, samples of oil being drawn from 
the transformers every few months for testing. Dielectric constant measure- 
ments were also made on these samples. 
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A sample of oil No. 2, oxidized in the laboratory, was used in the first 
series of power factor measurements. This oil solidifies between —50°C and 
— 60°C. Another oil which solidifies at about — 10°C was oxidized and used in 
further experiments. In order to increase its viscosity, it was mixed with an 
equal quantity of a pure hydrocarbon lubricating oil of viscosity at 40°C, 18 
times greater and of about the same freezing point. 


RESULTS 


The data are all shown graphically. In Figs. 1, 2 and 3, the increase in 
dielectric constant of the oils is plotted against the time of heating. In Figs. 1 
and 2, curve No. 1 represents oil No. 1, curve No. 2 oil No. 2, ete. Similar 
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Fig. 2. Oils heated in atmosphere of nitrogen. 


data for the same oils, but obtained from measurements on samples with- 
drawn from transformers in service, are shown in Fig. 4. The data obtained in 
one of the chemical tests on these samples are shown for comparison in Fig. 
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Fig. 3. Oils heated in air in presence of copper. No. 1:—-Oil No. 1 heated with limited access to 
air. No. 2:—Oil No. 1 heated with free access to air. 


5. Figs. 6, 7 and 8 show the variation of dielectric constant and power factor 
with temperature for these oils. The unused oil of Fig. 6 is very similar to oil 


No. 1. 
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DiIsCUSSION OF RESULTS 


It may be pointed out that the dielectric constants of these oils when new 
give some indication of their composition. Aromatic hydrocarbons are known 
to have larger dielectric constants than naphthenic hydrocarbons, which in 
turn are larger than those of the paraffin hydrocarbons. For instance, the 
dielectric constant of benzene, an aromatic hydrocarbon, at 30°C is 2.26,; of 
cyclohexane, the corresponding naphthenic hydrocarbon, 2.013, and of hex- 
ane, the corresponding paraffin hydrocarbon, 1.87;. The dielectric constants 
of the new oils are 2.241, 2.159 and 2.110 for No. 1, 2 and 3 respectively. Oil 
No. 1 containing an appreciable percentage of aromatic compounds, has the 
largest dielectric constant, and No. 3, the paraffin base oil, has the smallest 
value. 
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Fig. 4. Change of dielectric constant of oils in use in transformers. 


The oils show a marked difference in the development of polar character- 
istics. The final oxidation products under these conditions are chiefly acids 
and esters. The formation of a certain amount of solid matter or sludge also 
occurs. The oxidation products are all polar and have relatively large dielec- 
tric constants. Of course, as these compounds are formed in the oil, the dielec- 
tric constant of the mixture increases. That the presence of the oxygen in the 
air was the most important factor in the destruction of transformer oils in 
service is well known, and is brought out clearly here by the striking differ- 
ence in the development of polarity in the oils when they are heated in air and 
when under an atmosphere of nitrogen, as shown in Figs. 1 and 2. The cataly- 
tic influence of copper is very clearly shown in Fig. 3. Oil No. 1 was used in 
this test and hence curve No. 2 is comparable with curve No. 1 of Fig. 1. [tis 
readily apparent that ‘the presence of the copper almost doubles the rate of 
decomposition. The change in the dielectric constants of the oils as they de- 
velop polar characteristics in service is shown in Fig. 4. These curves are very 
similar to those of Fig. 1, and indicate that the results obtained in the lab- 
oratory tests closely parallel those obtained in service. The polar substances 
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formed by heating oil exposed to air are practically all oxidation products, 
chiefly acids and esters. A chemical test which measures the amount of acids 
and esters formed should show at least rough parallelism to the change in 
dielectric constant. Such data are shown in Fig. 5 and the similarity to Figs. 1 
and 4 is apparent. A further proof of the polar nature of the used oils is shown 


SAPONIFICATION VALUE (Mg KOH /GRAM OF O1L) 
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Fig. 5. Development of spaonifiable matter in transformer oils with service. 


in Fig. 6. Here the dielectric constants of the oils are shown plotted against 
temperature. Their dielectric constants do not differ much at the extremely 
low temperatures where the oils are solid or nearly so, but the used oils show a 
definite maximum at somewhat higher temperatures where the oils become 
fluid again. The unused and supposedly nonpolar oil shows a very slight de- 
crease in dielectric constant as the temperature is lowered below — 50°C. This 
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Fig. 6. Variation of dielectric constant of transformer oils with temperature. 


may be due to the contraction which occurs when the oil soldifies, with the re- 
sult that the space between the condenser plates is not completely filled. 

The variation of power factor with temperature as shown in Figs. 7 and 8 
indicates that it is very largely caused by dipole energy loss as the molecules 
rotate with the field. The relaxation time, that is, the time in which the or- 
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derly arrangement of the molecules caused by an electric field will be reduced 
to 1/e of its initial value is defined‘ by the equation ¢=82nr?/2kT where k is 
the Boltzmann constant, T the absolute temperature, 7 the radius of the 
molecule and » the inner friction constant of the liquid. By using the ordinary 
viscosity for the inner friction constant, it is possible to calculate a character- 
istic frequency at any temperature for a liquid of known molecular diameter. 


NO | 500000 CYCLES 
NO 2 1,000,000 CYCLES 


PER CENT POWER FACTOR 





- + 
TEMPERATURE - °C 


Fig. 7. Variation of power factor of oxidized transformer oil with temperature and frequency. 


At frequencies either very much above the characteristic or very much be- 
low it, the dipole contribution to the power factor is obviously quite small. In 
the present case, instead of measuring the power factor at constant tempera- 
ture over a wide range of frequencies, the more convenient method of varying 
the temperature was used. Starting at about 100°C, the power factor is practi- 
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Fig. 8. Variation of power factor of oxidized transformer oils with temperature and viscosity. 


cally independent of temperature down to 20°, where a sharp increase begins. 
As the temperature is lowered still further and the viscosity increases rap- 
idly, the power factor passes through a maximum and drops back to a very 
small value. When the dipoles are so firmly held that they can no longer fol- 


' P. Debye. Polar Molecules, p. 84. Chemical Catalog Company, 1929. 
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low the field, their contribution to the power factor of course becomes negligi- 
ble. The higher the frequency, the higher the temperatue at which the maxi- 
mum should occur, and we find the maximum of the 1,000,000 cycle curve 
about 10° higher than that of the 500,000 cycle curve. 

The most interesting feature of Fig. 8 is that in curve No. 1 the maximum 
occurs at —40°C, about 30° below the freezing temperature of the oil. This 
shows that although the oil is quite solid, the molecules are still free to rotate 
with the field, and hence that the mechanical viscosity and the inner friction 
which opposes the turning motion of the molecules are not always the same. 
The peak in curve No. 2 is very much broadened and it is shifted from — 40° 
to —10°. The power factor of this more viscous oil is less, due to the dilution 
with a non-polar oil, and the maximum is shifted to a higher temperature be- 
cause of the increase in viscosity. Johnstone and Williams® have shown that 
molecular radii of the right order of magnitude are obtained when the ex- 
perimental values of viscosity and characteristic frequency are substituted in 
Debye’s formula for the relaxation time. However, the fact that the maxi- 
mum in power factor (and in dielectric constant, although the data are omit- 
ted for brevity) occurs many degrees below the solidification temperature of 
certain oils shows that the inner friction constant or molecular viscosity is not 
always even approximately equal to the mechanical viscosity. 

The measuring of the polarity of transformer oils in this fashion should 
serve as a useful indication of the extent of decomposition in a given used oil 
and as a convenient test by which to predict the behavior of oils in trans- 
formers. It gives no indication of the tendency of an oil to form sludge except 
insofar as the degree of oxidation is an indication of the tendency to sludge. 
The presence of suspended sludge in the oil does not appreciably affect its 
dielectric constant. However, the method does offer a quick and compara- 
tively easy means of differentiating between several oils as to their tendency 
to decompose. 

The author is pleased to acknowledge his indebtedness to Mr. J. G. Ford 
for the oils studied and for the data in Fig. 5, and to Dr. C. F. Hill for his 
helpful advice and criticism throughout the course of the work. 


* Johnstone and Williams, Phys. Rev. 34, 1483 (1929). 
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Potential Distributions in Large Cylindrical Disks with 
Partially Penetrating Electrodes 


By Morris MusKat 


Gulf Research Laboratory 
(Received March 3, 1932) 


Because of the equivalence of the pressure in a liquid bearing sand to a velocity 
potential, an analysis has been carried through of the potential theory of the electrical 
analogue of an artesian well partially penetrating a water bearing sand. 

Part I. In Part I the method of analysis is illustrated by the solution of the prob- 
lem of a hemispherical electrode imbedded in one of the faces of the disk. Two types 
of formulae are developed suitable for study of the potential distribution at small and 
large distances from the electrode respectively. The resistance of this system as a func- 
tion of the disk thickness is computed and plotted. 

Part II. In Part II is treated the more general problem of the potential distribu- 
tion for a partially penetrating electrode, with a uniform flux density along the elec- 
trode surface. Separate formulae for the potential function are again derived that are 
convenient for discussion at small and large separations from the electrode. It is shown 
that they reduce in the limit of small penetrations to those of Part I and in the limit of 
complete penetration to well-known radial flow formulae. Expressions are also given 
for the stream functions corresponding to the derived potential functions. Finally it 
is shown that when the flux density along the electrode is assumed to be constant the 
potential over the electrode is far from uniform, falling sharply near the electrode ex- 
tremity to a value slightly more than half of that at the top of the electrode. 

Part III. The problem of physical interest in which the electrode is at uniform 
potential is attacked in Part IIT. After discussing the case of a flux density over the 
electrode surface increasing with distance from the top the direct method of superpos- 
ing discontinuous flux density elements is finally used to obtain electrode surface” 
potentials uniform to within 2 percent of the average surface value. The equivalence 
of this method to the direct solution of the integral equation for the problem is pointed 
out. The resistances of the disk for various electrode penetrations and disc thicknesses 
are then computed and plotted. From these a semi-empirical approximation is derived 
permitting the resistances for other cases to be derived from a single expression with- 
out the necessity of first adjusting the electrode flux distribution so that its potential 
be uniform. 


LTHOUGH the problem of the potential distributions in cylindrical 
disks with partially penetrating electrodes appears to have no practi- 
cal interest electrically, it is of considerable importance when applied to the 
question of flow of liquids through sands into artesian wells. This significance 
is due to the fact that the fluid pressure in a sand through which a liquid is 
flowing, being equivalent to a velocity potential, obeys Laplace’s equation. 
To find the pressure distribution in a sand one has therefore to solve the 
corresponding potential problem for the same geometrical system. 
This analogy may be formulated more precisely in the following way. For 
a sand bearing a liquid in motion, the velocity potential, ®, is related to the 
pressure Pp, at the point (7, 0, z) by the equation 


& = k(+ p — 182), (1) 
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so that 


OP 1 of O® 

Sa -79;% = -—--——! &2e-—- =: & e- =») (la) 
or r 06 Oz 

v° = 0. (1b) 


In these equations, k is the permeability constant of the sand for the 
fluid in question—of density y—, and is defined best by Eqs. (1) and (1a) 
themselves; g is the acceleration of gravity and 6 is the “effective” fluid 
velocity vector, interpreted as the volumetric fluid flux across a unit section 
of sand, its components in cylindrical coordinates being v,, ve, and v., ve in 
the present discussion, always vanishing by symmetry. 

Eqs. (1), (1a), and (ib) are really nothing more than the analytical 
statement of D’Arcy’s law established experimentally as an empirical sub- 
stitute for the Stokes-Navier equations, whose application to the flow in a 
sand would be a hopelessly difficult and complicated one. Eqs. (1a) further 
state that v, and v, are positive when directed towards positive 7 and 2, the 
latter being taken as downwards. 

It is clear that electrically the sand corresponds to a resisting medium 
of specific conductance k, voltage distribution , and current flux distribution 
d. Although, from the electrical point of view there is no analogue for the 
gravity term ygz, still, since the problem is sufficiently defined by the poten- 
tial function ®, there is no need to specify—except in the final interpreta- 
tion—whether it may be split up as in Eq. (1) or not. The electrical analogy 
also implies that a well penetrating a fluid bearing sand corresponds to a con- 
ducting electrode and the impermeability of the sand faces corresponds to 
electrical insulation. Since the present discussion will be essentially a poten- 
tial theory analysis of the problem in question, the electrical phraseology will 
be used throughout. 
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Fig. 1. Diagrammatic representation of disk of radius 7; and thickness a/2, penetrated by an 
electrode of radius rz and length b. 


STATEMENT OF THE PROBLEM 


. 


For simplicity the analysis will be confined to the case of cylindrical disks 
faced by insulating parallel planes, i.e., such that over them the normal 
gradient of the potential vanishes. At the center of the disk is inserted a con- 
ducting electrode of small radius rz, and length b. Corresponding to the case 
of practical interest when applied to artesian wells, it will be supposed that 











POTENTIAL DISTRIBUTIONS IN CYLINDRICAL DISKS 331 


the outer radius of the disk, 71, is very large as compared to the electrode 
radius, ’2, and of the same order as or larger than the thickness of the disk, 
a/2 (cf. Fig. 1). 
Analytically stated, the problem is: To find a solution of the equation 
eb Ph Ph 


v2) = + — + —. — 0, (1b) 
Ax? © Ay? © dz? 





such that: 


O® 
-=0,2 = 0, a/2; (1c) 
Oz 
O® 
— = const. 
or r=f2,2<); (1d) 
or: ® = const. 
® = const; r = rj. (le) 


Condition (le) will not be strictly adhered to. Since, as was stated, r; 
will be taken as much larger than 72 it follows from elementary considerations 
that necessarily ® will be very approximately constant at r=7;. For all cases 
of practical interest the slight variation of ® with z at r=r, will be of no im- 
portance. Hence in the following analysis no attempt will be made to force 
® to have a strictly constant value at r=r, but rather, the ® that satisfies 
the other conditions will be taken as satisfying with sufficient accuracy the 
last condition also. 

It should be mentioned, too, that in the preliminary development of the 
analysis, as is given in Part II, the first of conditions (ld) will be satisfied. 
That is, a solution will be developed and discussed, for which the flux density 
along the electrode is uniform; from this, a solution will be developed in Part 
III such as will satisfy the second of Eqs. (1d) to a close approximation. It is 
true that only the latter solution is of physical interest, but the first case will 
be solved essentially as an analytical preliminary problem. 

From the symmetry of the problem it is clear that one should use cylindri- 
cal coordinates. It will be convenient to place the origin at the center of the 
upper face and to take z positive as directed downward, as already men- 


tioned. It will also be convenient to introduce from the outset the dimension- 
less variables: 


= W;— = X; — = p23} — = pi. (2) 


Before entering on the analysis it should be mentioned that the present 
problem has recently been investigated by A. F. Samsioe.' Since he was not 
interested in details of the potential distributions, he derived formulae which 
are not suitable for studying this phase of the problem. Thus, for example, 
at r=0, his expressions appear to diverge even if z>b. Further, for large 
values of r his equations are rather awkward to handle. Finally, Samsioe has 
not attempted to satisfy in detail the requirement that the electrode be at a 
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uniform potential, but treated only the case where the flux density along 
the electrode is uniform. 

Because it has some features not characteristic of the general penetrat- 
ing electrode and because it turns out to be the most important case for an 
application made in another paper, it will be useful to discuss first in detail 
the special problem where the electrode does not penetrate the disk but is 
embedded in its upper face. This problem will therefore be solved first. 


I. THE NONPENETRATING ELECTRODE 
A. Potential function for small values of p 


By a nonpenetrating electrode is simply meant a semispherical cap em- 
bedded in the upper face of the disk. For purposes of the analysis, however, 
it will be convenient to replace the electrode by a small spherical equipotential 
surface enveloping a mathematical source of strength g, placed at the origin. 
The fundamental solution of the potential equation will therefore have the 
form: 


0 ee : (3) 


(r? + 3?)1/2 a(p? + w)!/? 








To obtain a solution satisfying Eq. (1c) the well-known method of images, 
also used by Samsioe, will be employed. Thus, it will be readily verified that 
if images of the source at the origin be placed along the axis at the points 
(0, tna) the resulting potential distribution will give zero flux across the 
planes z=0, a/2. With this arrangement the problem is formally solved, and 
there remains only the development of suitable formulae to represent the 
solution. 

Thus, from Eq. (3) it is seen that the complete solution has the form: 

- x0 
= £ an — easton —— + 2. yon + const. (4) 
a 0 (p> + (n+ w)?)!”? a (p? + (n — w)2)!2 





Looking at Eq. (4), without reference to Eqs. (1c) and (1d), it will be noted 
that it represents the potential due to an infinite regular array of equal point 
sources or charges, lying on the z-axis, and spaced by the distance a. From 
this point of view Eq. (4) has been discussed by Madelung? as a preliminary 
to his analysis of the potential fields in space lattices of heteropolar crystals. 
Although Madelung does not concern himself, in the case of the line array, 
with details of the potential distribution, it will be convenient to make use of 
some of his transformations. 

As is to be expected of the potential due to an infinite set of like images, 
Eq. (4), as it stands, is everywhere divergent and appears therefore to have 
no meaning. However, this difficulty is only apparent and of no physical 
consequence since one may subtract an arbitrarily large constant from the 
summations and in the limit obtain a convergent representation of the 
variable field at all points of interest. This will appear more clearly presently. 


1 Samsioe, Zeits. f. ang. Math. u. Mech. 11, 124 (1931). 
2 Madelung, Phys. Zeits. 19, 524 (1918). 
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Thus, neglecting for the moment the question of convergence, the series 
in Eq. (4) can be expanded for: p<1—w21/2, i.e., for small values of r, as: 





! I< I 1 p* 
@ = —— pee + — yon 7 = yy —_ 
a(p? + w?)!/? a n+ w \ 2 (n+ w)? 
3 p* , 2 1 1 p | 
~~ ——_—— + - yon yy _ : (5) 
8 (n+ w)! a 1 n— wl 2 (n— w)? | 
3 p* ) 
es wee ss t + const. | 
8 (n — w)! 
Now, from the theory of gamma-functions one has the relations: 
I’(y) 1 . 1 1 
v(y) = —— = —6 —-— + lim Dm (—- : ) 
r'(y) y n—2 4 m yt+m ; 
where , (6) 
6 = lim ( yom — — log ") = (0.5772. | 
"0 i m ) 
Hence: 
: 1 ‘ 1 a: tea 
[f- <teeasenne 2 =- oe = - im log 7, 
1 n+ w w n> 2 | 
(6a) 
rs 1 | 
n ——— = — (1 — w) + lim logn. )} 
1 = | oe 
Further, by the definition of the ¢-functions, one has: 
— 1 1 
n ———— = ¢((s,w) - ; 
, (n + w)* Ww 
f 6b 
io a) 1 | 
5 cnn we Se wh 
1 (n— w)* 
Eq. (5) can therefore be rewritten in the form: 
q 1 I ‘ 
» = (uu) - —y(1— w) - 072 (3, % 
a L(p? + w’)!/? W 2 { 
1 \ 3 j 1 _ 
~— + (3,1 — w)$ +— 92505, w) - - (7) 
w* f gs | ti 
+ ¢(5,1 — x)! —---+2 lim log | + const. 


From this expression it is seen that the divergence of the series in Eq. (4) 
may be separated out into the term 2 lim, ..,. (log ») which can be formally 
removed by making the arbitrary constant term equal and opposite to it. 
It will therefore be dropped hereafter although a finite constant term may 
be added so as to make possible an arbitrary choice of ® at some given point. 
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Since the y-functions are usually tabulated only for arguments greater 
than 1 it is useful to transform the functions in Eq. (7) so that they satisfy 
this condition. This may be done by means of the well-known relations: 


1 
viv) = —-—+¥7(14+ y); VU — vy) = vy) + oreetn sy. (6c) 
Bs 


Applying these to Eq. (7) one gets finally: 











q . 1 1 ) 
&®=—)| — Wit w)+— — retnew+ : 

a | Ww (p? — w?)! 2 
Bf 1 

ae iat w) — co + (3, i- WwW) f (8) 
3 f 7” 1 " 

+ —- pf ¢(5, yo-- + ¢(9, 7 w) oF ° 
8 \ ws : 


It is readily shown that for very small distances from the origin Eq. (8) 
reduces to: 





> (p? + w?)'? €1. (8a) 


This means physically that very near the origin the potential is spherical 
and due almost entirely to the single source at the origin, and this is, of 
course, what is to be required if the electrode is to have a spherical shape. 

From Eq. (8a) one can readily obtain the total flux Q leaving the electrode 
and flowing through the disk. Thus: 


a’ 
Q= f- do = + 2rq, (9) 


on 


where in the integral, 7 is the normal at the element do on any surface en- 
closing the origin, the integral extending, however, only over that part of 
the surface lying between the planes z =0, and z=a/2. Eq. (9) simply verifies 
our interpretation of gas the source “strength”. 


B. Potential function for large values of p 


For small values of p Eq. (8) can be handled numerically without diffi- 
culty. For large values of p, Eq. (8) diverges; but even the corrected series— 
obtained by separating out terms like 1/(p?+w*)'? from the series in Eq. 
(4)—is quite unsuitable for detailed study. In fact, one must abandon the 
power series in p. As a substitute it will be convenient to make use of Made- 
lung’s rapidly convergent series obtained by a direct solution of the potential 
equation for a continuously distributed charge along the z-axis, and so ad- 
justed that it goes over into the potential due to a point array on properly 
passing to the limit of a concentrated charge density. 











POTENTIAL DISTRIBUTIONS IN CYLINDRICAL DISKS 335 


The equivalent of Eq. (8) becomes, by Madelung’s representation: 


_ ) 
p = “\4 >on K o(24np) cos 2rnw + 2 log | (10) 
a l p 
where Ky is the Hankel function of order 0. For p>0.8 the first term in the 
sum suffices for almost all purposes. Again, in Eq. (10) an infinite constant 

has been dropped as it was in Eq. (8). 

It is readily seen from Eq. (10) that ® becomes practically independent 
of z as soon as r becomes of the order of a, i.e., the flux is essentially radial 
except in the region near the electrode no more distant than about twice the 
thickness of the disk. This will be seen more clearly from a plot of ® to be pre- 
sented later, and it justifies in detail the remarks made on page 331, that 
solutions satisfying Eqs. (1c) and (1d) will necessarily satisfy Eq. (1e) with 
sufficient accuracy. 

It is of interest to derive here the relation between the total flux Q and 
the potential difference A® existing between the electrode and some large 
distance r; from it. This may be found as follows. As: 


. 7 1/2 
Koy) 2(—) ev 
2y 


for large values of y one may omit the summation in Eq. (10) for large p; 
whence it reduces to: 
2q 
& = — log 2/p. (10a) 
a 
Now it has already been noted that for small distances from the origin ® 
depends only on (p?+w?)/?. Hence specifying that the electrode radius is 
given by (p?+w”)!/?=pe, the potential of the electrode may be expressed by: 


2 
© = O(w = po, p = 0) = 4| - 2V(1 + po) + — — retn ros]. 
a pe 
The potential drop then becomes: 
—¢ 2 2 
se - —42 log — + 2¥(1 + px) — — + wctn apa]. (11) 
a Pi p2 


Now A®/Q is clearly the electrical resistance of the disk for unit con- 


ductance of the disk material. Denoting this resistance by Q, it follows from 
Eqs. (11) and (9) that: 
A® 1 2 2 
Q = — = — | — — rctn ape — 2f(1 + pe) — 2 log— }. (12) 
Q 27a L peo Pi 


Noting that for peX1, 7 ctn mp2, =1/p2 Eq. (12) may be approximated by: 





1 1 
R= — —(y¥(1) + log 2/p:). (12a) 


2tre wa 
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Hence, since the last two terms are evidently smaller than the first, ex- 
cept for small values of a, the resistance for the nonpenetrating electrode is 
practically independent of the thickness of the disk. To show this more clearly 
the value of 2 is plotted in Fig. 2 as a function of a for the case where r2 = 0.01 
and r,;=10. The range of a is from 1/2 to 10 so that pe. varies from 0.02 to 
0.001. 
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Fig. 2. Resistance of disk with a nonpenetrating electrode as a function of the thickness 
of the disk. Electrode radius r2=0.01 and disk radius r; = 10.0. Horizontal dotted line represents 
resistance of a semi-spherical shell of inner and outer radii of 0.01 and 10.0 respectively. 


One may note, too, that the last terms in (12a) simply represent the de- 
viation of Q from that for a perfectly semispherical resisting shell, since the 
resistance of the latter may be easily shown to be: 


re 1 1 
2 ie eel ee ae —) a (12b) 
2r \ re ry 2tre 


The horizontal line in Fig. 2 is the value of Q2 for purely semispherical 
flow and again with r2=0.01, r;=10. It is seen from this that as long as p» 
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Fig. 3. Equipotential curves in a disk with a nonpenetrating electrode placed at center of 
its top surface. 4, is the potential function of equation (13). p=radial distance /2 (disk thick- 
ness). w = (depth in disk) /2 (disk thickness). 


=re a<0.005 the resistance of the disk differs but slightly from that of the 
corresponding hemisphere. 
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To see more clearly the nature of the potential distribution both for 
small distances from the electrode and distances at which the radial char- 
acter of the flux begins to predominate, a plot of several equipotential curves 
in a plane through the z-axis is given in Fig. 3. To make the distribution in- 
dependent of any particular thickness of the disk the latter was eliminated 
by taking g/a=1 in Eq. (8). The particular electrode radius was also elimi- 
nated by leaving the absolute potential over the electrode quite arbitrary. 
Thus the potential function plotted in Fig. 3 is given by 


1 1 
&, = ———— + — — rctnrw — W1+ w) 
(ep? + ww)? w 
1 , i rer 
— riatl ; w) — a ¢(3, 1 —_ w) ¢ —_ 8 p (5 uw) 
(13 
1 





i ) 
®, = 4 > nK o(2enp) cos 2rnw + 2 log 2/p:p > 3. 
1 
Although tables of the ¢-functions were available for use in the first value 
of ®; for p<1/2, it was found just as convenient to use the second form with 
three terms in the sum even for p = 0.3. 
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Fig. 4. Diagrammatic scheme of image system used to represent potential functions for par- 


tially penetrating electrodes. 


It will be noted that, as already indicated before, although near the elec- 
trode the equipotential surfaces are spherical, just as about a single point 
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charge, 2s soon as one recedes from the electrode to a distance less even than 
twice the thickness of the disk the equipotential surfaces become closely 
cylindrical as if the electrode penetrated the disk completely. 

For the present purposes the above discussion will suffice for describing 
the features peculiar to the special problem of the nonpenetrating electrode. 
As to the other details, it will be better to bring them out in comparison with 
the cases of penetrating electrodes. These will therefore be analyzed in the 
following section. 


II. PARTIALLY PENETRATING ELECTRODES WITH 
UNIFORM FLUX DENSITY 


A. Potential function for small values of p 


To analyze the case of the partially penetrating electrodes, it is convenient 
to begin with the potential function due to an element of source distribution 
located along the z-axis at a distance A from the origin. The elementary 
potential function in this case is of the form: 


gd qdu 
(P+ (@— dH) E+ (ww) 


where gd) is the strength of the element of length dX, and yu is the dimension- 
less variable \/a. 

Again, to obtain a function ® which will satisfy Eq. (1c) it is necessary 
to add to the original source at (0, \) an appropriate set of images. These, 
it will be readily verified, should be placed at the points (0,.+n”a+d). One 
therefore has, for the differential potential function due to the element at A 
and its images, the value: 








(14) 


1 1 
: ae aos ; 
(p? + (w— m)*)'? — (p? + (w + w)*) 


at 1 ‘. 1 
1 (+ (n+ w—p)?)? (0? + (n+ w+ p)?)'? 








d& = adn 
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1 1 
+— — | } |. 
(p? + (2 — w— p)*)"* = (p? + (2 — w+ p)*)"?) 


* To be exact the first term in Eq. (15) should be changed to 1/(p?+(u—w)?)'? for w< uy. 
However, this is of no significance except when p=0. In the following, therefore, Eq. (15) will 
be used as it is, since a detailed check will show that all the resulting formulae are rigorously cor- 
rect. It should be pointed out, too, that the author does not pretend that the formal solution 
in terms of images, as expressed by Eq. (15) is new or original for as early as 1873 Maxwell set 
up this same system for a point electrode in a disk and even for the more general case where 
the disk is surrounded by other conducting media. Further, geophysicists interested in electrical 
prospecting have recently worked on this more general problem of Maxwell. However, no de- 
tailed development and discussion of Eq. (15) or of its equivalent for large values of r as pre- 
sented here has been found in the literature, except in the references specifically cited, nor has a 
treatment of the partially penetrating electrode been given by anyone except by Samsioe, as 
already mentioned. 
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Supposing now that the source at A is continued with constant strength 
along the z-axis from y=0 to y=8, i.e., from uw =0 to p=x, the potential due 
to the continuous and uniform distribution will be given by: 








b 
o= [as 
0 
- af , a np oat w+ xt (p? + (n — w+ x)?) ise (96) 
; n—-w—xt(e+(n—w— x2 | ; 
~ ( ( w : 
» Die oe staid Beth a S| 
r n+w—x+(p?+ (n+ w-— x)?)! ) 


This is essentially the procedure used by Samsioe. However, his counting 
of the images led him to log terms whose arguments vanish for p=0, al- 
though one could easily transform his expression so as to be equivalent to 
those in (16), before passing to the limit p =0. In Eq. (16), ® has a singularity 
at p=0 only for w<x, for then the term log {w—x+(p?+(w—x)?)"*} 
hence ® becomes infinite as log p, as indeed it should. 

Nevertheless, Eq. (16) is not particularly convenient to handle unless 
one goes through lengthy expansions of the individual logarithmic terms. 
It is therefore better to expand the series in Eq. (15) directly, as was done 
in the case of the nonpenetrating electrode and then integrate. By this 
method, then, one obtains first: for p< 1—x—w, to terms in p*, 


and 
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The restriction p<1—x—w can be advanced to p<1—w or p<1—x by 
removing in the first summation in Eq. (17) the terms for = 1 and separating 
them out into the single term 1/(p?+(1—w—y)?)"”, but since in the applica- 
tion and discussion to be made later the series expansion for ® will be used 
only for p<1, it will be somewhat easier to use Eq. (17) as it stands. 

Now by the aid of Eqs. (6a) and (6b) for the y and ¢-functions, Eq. (17) 
may be written as: 
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Again, the divergence of the series in Eq. (15) appears to be in the term 
4 lim,.. log m, and for the same reasons as already discussed this term will 
be dropped in the following. 

As the potential d® will have to be integrated with respect to u the various 
terms must be expressed in differential form. This may be done on noting the 
following relations: 


YVi-—-w-—p) =VW(wtyu)tretnarwt yp), 
vViw— vw) = WI — wtp) — retna(w — p), 


d | (18) 
v(m + vw) = — log I'(m + wp), 
du | 
7 1 0 1a) 
¢(s, w).= — —f(s— 1,4). 
1—s Op 


Integrating now Eq. (17a) and making use of the well-known relation: 
us 
r(w)l(l — wv) = ——» (18a) 
sin rw 
one obtains finally: 


= Oo ) 
f — dy 
0 Ou 


af 2 log r — 2log l'(w+ a)F(1 — wt x) 
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— log sin r(w + x) sinx(w — x + log ——— 
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This is the final form for the potential; it is useful numerically for small 
values of p and valid only for constant flux density, q, along the electrode. \t will 
be discussed in detail in sections C and D. For large values of p another ex- 
pression will be derived presently. 


waxt(eet+(wtx))? 1 oy 
2)! 





Je 





B. Potential functions for large values of p 


Since Eq. (19) is essentially a power series in p, it is clear that it is not 
suitable for studying ® for large values of p, i.e., at large distances from the 
electrode. A suitable expression can be obtained, however, by the following 
extension of Madelung’s expression for the potential due to a point lattice of 
equal charges, as is given in Eq. (10). Noting that the system of images as 
indicated in Fig. 4 is equivalent to a superposition of two such lattices, it is 
readily verified, on comparing with Eq. (10), that the potential due to the 
source at \ and its set of images, is given by: 


db = gqdu E yon Ko(2rnp) } cos 2xn(w — p) 
1 
* 
_ (20*) 
+ cos 2rn(w + u)} + 4 log - | 
p 


— 


* The extended character of this whole analysis may seem artificial when it is observed 
that a direct solution of Laplace’s equation for the present problem gives the simple results for 
a flux element gdu at u: 


ws M 
7 * Jo(pd) cosh Aw cosh A(4 — pw 
> = 2au)du f nn te ment = — 3-4), 
. sinh — 
2 
wu 


* Jo(pd) cosh Au cosh A(4 — w) 


ae a 
sinh — 
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@ = 2gluddu f 


However, a little inspection will show that a detailed discussion of these integrals would 
lead to an analysis and results quite equivalent to the lengthy developments used here. In fact 
it is not difficult to show by contour integration that the above integrals have an expansion in 
terms of the Hankel functions identical with Eq. (20). The use of these integrals as a starting 
point for a detailed discussion of the problem would therefore not shorten the work over the 
somewhat more direct and more familiar method of images used in the text of this paper, for 
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Integrating now from 4=0 to w=x, with the assumption that g is con- 
stant, one readily finds as the equivalent of Eq. (19): 


; & Ris P 2 
® = 4g| — >on — Ko(2nnp) cos 2rnw sin 2xnx + xlog—|. (21) 
T 1 n p 
Again, as the Hankel function Koy decreases exponentially for large argu- 
ments, Eq. (21) is very easy to handle numerically even for relatively small 
values of p. Eq. (21) also permits an easy calculation of the total flux through 
the disk, which is given by 


19 


‘ a® ' f° ae sili 

= —do = — /zé o—— 6D, 2 
on 0 Op 

where 7 is the normal to the element do of the surface enveloping the elec- 

trode. The second integral implies that the enveloping surface is a cylinder 

coaxial with the electrode. Applying then Eq. (21) to (22) it is seen that the 

trigonometric terms give no contribution and the result is obtained that: 


Q = 4raqx = 4rqb, (22a) 


which again is really nothing more than a confirmation of our physical inter- 
pretation of the coefficient g. 

Eq. (21) shows further that for appreciable distances from the electrode 
the potential distribution, just as in the case of the nonpenetrating electrode, 
takes on a radial character and loses its dependence on z. And as is to be ex- 
pected, this approach to a radial flow character is more rapid in the case of 
the penetrating electrode. This is seen from Eq. (21) in the fact that each 
term depending on w, and hence on 2, has a factor sin 2n7x/n <1 in addition 
to the Hankel function which alone occurs in the corresponding Eq. (10) for 
a nonpenetrating electrode. 


C. Limiting cases of zero and complete penetration 


Before considering the details of the potential distribution it will be of 
interest to examine two limiting cases with respect to the electrode penetra- 
tion. Thus, first, it will be shown that it reduces to the well-known radial 
flow potential, independent of w, when the electrode completely penetrates 
the disk, i.e., when b =a/2, x =0.5. 

This may be done as follows: Noting that: 


1 
f(s, y) a c(s, y — 1) = 
and: (18b) 


log{ w+ (p? + w?)'/2} = ng 20 + ~— a a —_ >< Sor ere 
we 32w 





even the image expansion of Eq. (15) can be also obtained directly from these integrals on ex- 
pressing csch (\/2) as a series of exponentials and then integrating term by term. The elegance 
of the integral representation of # is therefore only apparent, at least as far as a detailed dis- 
cussion of @ is concerned. 
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one readily finds that the series in p in Eq. (19), involving the ¢-functions, 
reduces to the value: 


2 log 2 + log (w + 3)(3 — w) 
— log {w+ 3+ (po? + (w+ 4)%)"?} [4 — wt (p? + (F— w)?)1/}. 


With regard to the remainder of Eq. (19) one makes use of Eq. (18a) and 
then after a slight manipulation it reduces to: 


© = 2q log La sah. (23) 
p 
which is, in fact, the exact value that one should have for pure radial flow 
from a uniform continuous source along the axis of the disk. 

The other limiting case of interest is that where the penetration of the 
electrode is very small and approaches the case of the nonpenetrating elec- 
trode. It will therefore be shown now that for x very small, the general Eq. 
(19) does reduce to the expression already found for the nonpenetrating 
electrode. To do this, it is useful to note first the following relations, valid 
for small values of x: 


log '(y + x) = log I'(y) + wy); 


log sin r(y + x) = logsinry + xamctn ry; 








ytxt (p? + (y+ x)?)'? ytu 2x 2x 
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t(s,y + x) = F(s, y) F sxt(s + 1, y). 


Putting these equations into Eq. (19) and again making use of Eq. (6) 
one obtains finally: 











1 1 
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The expression in the bracket is clearly identical with the corresponding 
expressions for the case of the nonpenetrating electrode as is given in Eqs. 
(7) and (8). Just as in this latter case the factor g/a represented the quantity 
Q/2zxa, where Q is the total flux from the electrode, so, here the factor 2qx 
is also equivalent to Q/ 27a, as is seen from Eq. (22a). 

The limiting values of the potential function are still more easily derived 
from the expression in Eq. (21) developed for use for large values of p. Thus 
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for the case of complete penetration one has only to set x =1/2 in Eq. (21) 
to obtain at once: 


& = 2¢ log 2/p, 


as was found from the series expression for ® (cf. Eq. (23).) 
On the other hand for the nonpenetrating case, i.e., x small, one need only 
note that 


sin 2umx = 2nmrx:x small. 


This may be assumed without appreciable error, even for m large since then 
Ko(2rnp) is anyway so small that the contribution of the whole term is 
negligible. With the above substitution Eq. (21) becomes: 


a ) 
@ = 2434 yon Ko(2urp) cos 2nrw + 2 log =|, 
p 
which again is identical with Eq. (10) taking into account Eq. (22a). 

Hence it is seen that in the limits of very small penetrations the general 
case of the partially penetrating electrode goes over continuously into that 
for the nonpenetrating case, and in the other limit of complete penetration, 
the limit of pure radial flow is also approached continuously. 


D. Stream line distributions 


It is often convenient, in potential problems, to study the stream lines 
directly rather than by inference from the equi-potential curves. It is there- 
fore of some interest to have available the explicit expressions for the stream 
functions corresponding to the potential function formulae, although they 
will not be discussed further in the present paper. 

In terms of the present notation the velocity components are related to 
the stream function ~(p, w) by the equations: 





1 dp 1 0d) 
Ur = ee : — Sooo — — 
a*p Ow a Op ; 
(25) 
1 oy 1 0 
° a’p Op a Ow 
so that, in terms of the potential function ®: 
ab + Ob | 
¥(p, w) = — ap -{ —— du =4¢ fox dp. (26) 
Op Ow 


The results of evaluating the integrals of Eq. (26) by means of previously 
developed formulae for ® are: Flux element gdyat (0, u); potential functions: 
Eqs. (17a), (20): 
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w— yp Ww 











V(p, w) = aad] —— | + O(p?) (27) 
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> 
\IV 


a) 
¥(p, w) = sagdu| w + 2p din Ki(2nmp) sin 2nrw cos duu | (27a) 
1 


Flux segment: g(u) =1 OS <x; potential functions: Eqs. (19), (21). 
p<i1 
v(x, p, w) = a[(p? + (w+ x)*)/? — (p? + (w — x)*)"/2] + O(p?) (28) 


p21 


II 


p 2. K,(2nmp) . 
V(x, p, w) = 4a | +— Son ——— sin 2nrw sin 2nrx (28a) 


T 1 nN 


Of course, as in the case of the potential functions, the trigonometric 
series for y, as in Eqs. (27a) and (28a), are strictly valid for all values of p, 
but only because they are not as convenient for numerical work when p<1, 
have the other expressions been developed. 


E. Potential distribution on electrode surface 


Up to the present the discussion has dealt with potential functions which 
by construction give no flux across the faces of the disk in which the electrode 
is imbedded and which have been shown to have very nearly uniform values 
over large cylindrical envelopes coaxial with the electrodes. In the derivation 
of these functions it has been explicitly assumed that the flux density, gq, 
along the electrode surfaces, was uniform. At the very beginning, however, 
it was pointed out that physically the electrode must correspond to an equi- 
potential surface, or that over its surface, ® was to be constant, as required 
by the second Eq. (id). The problem therefore remains of seeing to what 
extent the potential functions discussed above do satisfy this condition and 
how they may be made to do so in case they vary considerably over the 
electrode surface. 

In considering the potential distribution over the electrode surface, i.e., 
for values p=peo, w Sx, it is to be remembered that pe<1 for those electrodes 
that are considered of practical interest from the point of view of such appli- 
cations as the one to be discussed in another paper. With this in mind one 
can drop the series in p in Eq. (19). The singularities of the coefficients of this 
series for w=x are only apparent since: 


1 
(s, y) ne c(s, y + 1) 
ys 
= ¢(s,1) = 0(1) for y = 0;s 2 2. 


Since the whole series development of Eq. (19) is valid only if p<1—w—x 
the case w=x 1/2, for which ((s, 1-w—x) =(1—w-—x)~*, will cause no 
difficulty as these terms are cancelled by the equivalent approach to a singu- 
larity of log sin r(x +w) asx=w—1/2. 


' 
' 
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One can therefore safely reduce Eq. (19) to the form: 


pm = af 2 log r — 2 log T(w + x)I(1 — w+ x) 
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— log sin r(w + x) sin r(w — x) + log 
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Expanding the logarithmic terms and noting that 


1 pe” p 
log {| — b + (p22 + 6%)'/2} = “e-; : ars <1, 


one can express ® in the following way: 
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(29a) 





It therefore appears that in the limit of very small pz where the terms in 
Eq. (29a) independent of p2 are small compared to log pz, the potential at the 
bottom of the electrode is but half as large as that along the upper part of 


the electrode, i.e., 


2r 
@= 2glog— : x-—w>pe 
p2 
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2r 
@= glog—:x-—-we=0. 

pe 
Thus the potential along the electrode is far from constant. The variation of 
® along the electrode was noted by Samsioe from a numerical calculation 
for the case where the electrode had but one image—its own reflection in the 
plane z=0—but from the above it appears that in the limit the variation 
always involves a factor of 2. 

The physical reason for this factor of 2 is not difficult to see. Thus for 
points very close to the electrode, the potential at such points is largely due 
to the flux elements in their immediate neighborhood both above and below. 
At the end of the electrode, there are no more flux elements below the end 
point to contribute to the potential and its magnitude can therefore be only 
about half as large as that farther up the electrode where there are flux ele- 
ments on either side of the point to contribute to its potential. This will be 
shown in section A, Part III, to hold even for a variable flux density along the 
electrode surface. 
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Fig. 5. Potential distribution over surface p,.=1/600, as a function of depth, w, from the 
top of the disk; x = (length of electrode)/2 (disk thickness). Potential function # is that given by 
Eqs. (29a). 


When pz is not so small that the other terms in Eq. (29a) can be neglected, 
there is still a very considerable variation, although it does not reach the 
factor of 2. This is shown more clearly in Fig. 5, in which are plotted the 
functions ®/g, i.e., the potential for unit flux density, as a function of w for 
p2 = 1/600, for various values of electrode penetration x. Thus the curves up 
to w=x,—indicated by the horizontal lines—represent the potential dis- 
tribution on the electrode surface, the remainder up to x =0.5 being the dis- 
tribution in the disk just below the electrode, i.e., at the same distance from 
the axis. These curves show clearly that the assumption of a uniform flux 
distribution along the electrode surface must be dropped and an attempt to 
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obtain a uniform potential must be made by varying the flux density along 
the electrode. 


III. PARTIALLY PENETRATING ELECTRODES WITH UNIFORM 
POTENTIALS OVER ELECTRODE SURFACE 


A. Potential distribution on electrode surface for variable flux density 


It has just been seen that for a constant flux density, g, along the elec- 
trode surface the potential over the electrode is far from uniform. To elimi- 
nate this variation one should in principle consider g in Eq. (17a) or (20) asa 
function of u and determine it in such a way that after integration ® will be 
uniform for the preassigned value of p. This would lead to an integral equation 
of the form 


P *_ Ko(2nmp) ad ) 
—= }in ——— cos dnew f a(y) | 
4 1 nT 0 | 

(30)* 


2 . | 
sin 2mmudu + clog —> c = af pa(w)dp 
p 0 





in which a(u) is to be determined so that ® is constant for p=p2, wx. The 
exact solution of this equation would be quite difficult except for very small 
values of po, in which case the following treatment will indicate the nature 
of the result. 

For this purpose, it will be more convenient to use Eq. (20) directly. It 
will be assumed that g is everywhere finite and possesses a continuous deriva- 
tive. Then: 


p = if auddal >> Ko(2rnps) cos 2rn(w — p) 
0 1 





(31) 


2 | 
+ cos 2xn(w + w)} + log =| 
P2 


o= f q(u)dp. 


Then, as the above series is uniformly convergent as long as p2.>0, ® 
may be written as: 


Let: 


* It is of course appreciated that any function of the form: 


*cos 2nrwKo(2 1/2 
ou. Sala q(p2, #) cos 2nmudu + Blog p+C 
1 K o(2nxp2) z 


will give a strictly constant value for ® at p=p2 w <x. However, this value of & does not cor- 
respond to an electrode of depth x since the flux 6@/5p for p=0 does not vanish for w2x. One 
has therefore to choose q(p2, u) such that 6/4 =0 for w>x, and then the problem is essentially 
equivalent to the above in which the form of # ensures that 5@/5p =0 for w>x and the function 
a() is to be chosen so that ® is constant for w<x. 
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ps 2 2 1 
& = 40 log + > K o(24npe) E x) }sin 2an(x uw) 
p2 T 1 
+ sin 2rn(x + w } — | q'(u)dp} sin darn (u Ww) 


(3la) 


+ sin 2rn(u + w) |. 


Now Ko(p) has a logarithmic singularity at p=0. Hence considering ps» 
so small that log pz is the predominant term in the expansion for Ky it may be 
replaced by: 


Ko(2rnp2) = log (31b) 
pr 
even when 1 is large, since such terms anyway contribute but little because 
of the coefficient of 1/7 in Eq. (31a). ® may therefore be rewritten as: 


- 2 2 2 1 
® = 40 log + log > | g(x) }sin 2nr(x — w) 
p2 T pe n 
oa 
+ sin 2nar(x + w) } _ | q'(u)du}sin 2nw(u w) (31¢) 
0 
+ sin 2ur(u + w) |. 





Noting now that 
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_ 2 2 
© = 40 log — + 2(1 — 2x)q(x) log 
p2 pe 
2 2 1 . 
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and applying again Eq. (31d) one obtains finally,* 
2 
® = 2q(w) log — - (3 
p2 


2) 


For x =w, the sin 2mm7(x—w) terms in Eq. (31c) vanish. The result is then: 


2 = 1 7 1 . 
Q = 2 log — | 20 + g(x)(} — 2x) - —f q' (u) > —|{ — sin2rn(x— py) | 
p2 TT 0 nN (32 ’ 
‘ 32a 
+ sin 2nm(x + rau | = q(x) log — - 
p2 


It thus appears that for such small values of pe that the approximation 
of Eq. (31b) is justified the potential at a point on the electrode surface is 
simply proportional to the flux density at that point, the electrode extremity, 
however, requiring twice the flux density to give the same potential as that 
at a point on the electrode proper. Hence one may obtain even an arbitrary, 
distribution of potential on the electrode surface by simply assuming that the 
flux density g in Eq. (20) is proportional to the required potential distribu- 
tion, and doubling its value at the bottom of the electrode. 

However, as already pointed out, the solution just obtained of making the 
flux density distribution proportional to the required electrode potential 
distribution can only be valid if Eq. (31b) is a good approximation, i.e., for 
extremely small values of pz. In practical applications to artesian wells the 
well radius, 72 is usually about 1/2’ and the sand thickness a/2 is of the order 
of 50’, so that pp=1/200. Now log 2/p2=5.99 whereas Ko(2rnp2) <4.0 for 
even n=1. Hence it is clear that Eq. (31b) is not a valid approximation for 
cases of practical interest from the point of view of water flow into artesian 
wells. Before, however, presenting a method that will be applicable, it is of 
interest to see what a uniformly increasing flux density would give for the 
electrode potential. 


B. Potential distribution for a linearly increasing flux density 


Since the case of a uniform flux has already been discussed it will be 
perfectly general to take 
g(u) = mw. (33) 


Although the simplest formula would be obtained by applying Eq. (33) to 
(20), it will be better for numerical purposes to apply it to (17a) since it is 
particularly adapted for computation when pz: is very small. 

Denoting, then, the coefficient of g in the right side of Eq. (17a) by d® 
it is clear that for the present case: 


p~ = f ud®y = xPo(p, w, x) — f Po(p, w, udu, (34) 


0 


* If the same process is carried through for w>x it will be found that #=0. This simply 
means that g(u) =0, 1 >x as the form of (31) indicates. To get the true value of # for w>x, 
p2< <1 one must therefore include the constant terms in the expansion for Ko. 
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where by Eq. (19) 
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Evaluating the integral in Eq. (34) it is found that: 
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— (p? + (w — x)?)'/2? + 2(p? + w — (w+ x) log (w+ x) 


— (w — x) log (w — x) + 2w log w + ip? yw — w— x) 


Laas 
— Wi- w)+V¥1—-wtx)t+V¥(w — x) — Ww) | 
MMe +04 —— + — -=}4- “ot 43,1 | 
oa ; wt+e w—x w f Ve ren | 
— 2¢(3,1 — w) + ¢(3,1 — w+ x) + ¢(3, w — 2) | 
1 1 2 \ 
— 2¢(3, w) + §(3, w + x) — — — ee ae ome D 
(w+ x)? (w— x) wif 
+ I(p, w, x), 
where: 
— I(p, w, x) = f \2log F(w+ wl — w+ yp) | 
0 (34c) 
+ log sin r(w + yp) sin r(w — yp) \ du. 


Although it is difficult to integrate the individual terms in Eq. (34c) they 
may be paired and integrated in the following way. Thus, from the theory 
of gamma-functions it is known that: 
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2 log '(y) + log sin ry = log r + 2 Dona, sin 2nry;0< y< 1, (35)* 
1 


where: 
é, = 16 + log 2rn), (35a) 
un 
6 = 0.5772. 
Hence 


2) 


On ) 
—I=xlogr — > — cos 2ur(w + x) — cos 2urw} | 
1 “Tr 


> f {2 log P(1 — wt+uy) + log sin r(w — wu) j du. | 
0 


Since 1—w+y>1 for w>w, and (35) is valid only for y<1,—I must be 
evaluated separately for the cases of w>x and w<vx. For the latter case the 
integral in Eq. (35b) must be split up as: 


a ae 
0 0 w 


With this precaution one finds: 


w> x: 
x 
an { 
— I = 2x log r — 2 >) — cos 2nrw{ cos 2nrx — 1}; (35c) 
1 UT 
ws4 
a il 
n 
— I = 2x log r — 2 )> — cos 2nrw{cos 2nrx —1} 
1 NT 


f / (35d) 
+ 2(x — w)}log (x — w) — 1} + (x — w) log (— 1). 


It is to be noted that the term (x—w) log (—1) in (35d) is just right to cancel 
the equivalent imaginary term in (x—w) log (w—x) in Eq. (34b) for w<x. 
Hence ® remains real in all cases, as of course it should. 

By means of the above equations the potential distribution over the 
electrode surface, i.e., for pp =0.001 in the present case, was computed for a 
50 percent electrode penetration, i.e., for x =0.25. The result is indicated in 
curve I, Fig. 6. 

It will be noticed that the potential in curve I follows the flux density for 
about 80 percent of the length of the electrode and after that the end effect— 
in which the potential at the electrode extremity falls to about half of its 
normal value—begins to be significant. Strangely enough, however, when 
the proper uniform flux density is added to the linearly increasing flux one 


* It may be of interest to note that by using Eq. (35), after expanding Ko(2nzp) to terms 
in p?, one may show in detail the exact equivalence of Eq. (19) and (21) up to terms in p?. 
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obtains an almost exactly linear variation in the result potential. This is 
indicated in curve II of Fig. 6. 

In a manner similar to that used above other functional forms could be 
tried for g(u) and the corresponding potential distributions could be com- 
puted and plotted. It is clear, however, that one is not likely to find a suitable 














Fig. 6. Potential distributions over the surface of a 50 percent penetrating electrode of 
radius p2=0.001, for various types of flux distribution. Curve I: potential function of Eq. (34) 
due to a linearly increasing flux density. Curve II: potential function of Eqs. (34) and (29) due 
to a flux density g(u) = —1/4+ 4. Curve III: potential function (0.25, w) of Eq. (36) due to 
a uniform unit flux density. Curve IV: potential function A®(0.24, 0.25, w) of Eq. (36a) due 
to a unit flux density between w=0.24 and w=0.25; Curve V: potential function A&(0.15, 0.20, 
w); Curve VI: potential function #,(0.25, w) of Eq. (36b) due to a 0.01 unit source at w= 
0.25. 


form for g(u) unless complicated or “trick” functions are resorted to. It there- 
fore has appeared best to use the more direct, though less elegant, method of 
superposing flux elements at various places along the electrode axis and ad- 
justing their strengths so that the resultant potential over the electrode sur- 
face is approximately uniform. This method will now be described in more 
detail. 


C. Potential distribution due to a superposition of flux elements 


Since the discussion of this section will be confined essentially to the 
potential for very small values of ps—i.e., the potentials on the electrode sur- 
face—the terms in p,? will be definitely dropped from all the expressions and 
the potential functions used will be those given in Eqs. (29a). 

For definiteness it is convenient to introduce the following notation: 


2r } 
Oo(x, w) = 2 log — — 2 log I(x + w(x — w) | 
p2 Zz (36) 


— log sin r(x + w) sin r(x — w) } 
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A®(x1, X2, w) = Do(x1, w) — Pol xe, w). (36a) 
(x, w) = — Wl — w+ 2) — W1+ wt 2) 


—arctna(x + w) — retn a(x — w) 
1 1 1 


+> + + 36b 
xt+w x—w (po? + (w+ x)?)!/? —_ 


1 
(pa? + (x — w)2)"2 


Physically, the above equations have the following meaning. ®o(x, w) is 
the potential at (2, w) due to a uniform unit flux distribution on the z-axis 
and extending from 4» =0 to w=x. It is to be understood that the expression 
given for ® is to be changed as in Eq. (29a) if x—w~pe or x =w. AP(x1, x2, w) 
then represents the potential due to a unit flux element between w =x, and 
=xe. Finally, ®,(x, w) is the potential due to a unit source at up =x. It is ob- 
tained most conveniently from Eq. (17a) by setting gdu=1, w=x. Of course, 
all the functions ®), A®, and ®, satisfy Laplace’s equation and give zero flux 
over w=0, 1/2. 

The method of superposition is simply that of superposing various com- 
binations of the functions in Eq. (36) such that the resultant gives an ap- 
proximately uniform value for wSv~x. It is a rather tedious task and takes 
considerable time but after a few trials one’s experience will suggest proper 
“guesses” for the proper strengths and distributions of the functions. 

To indicate more clearly the nature of these adjustments and the mean- 
ing of the functions, plots of some representative cases, for pp=0.001, are 
given in Fig. 6. Although finally, the resultant potential function is to be ex- 
pressed as a sum of terms as %o(x, w) for different values of x, it is more con- 
venient at first to use the functions A®, since from their form one can more 
easily tell just how to build up the potential where it may deviate from its 
average value. The forms of the curves for A® as given in curves IV and V 
also show clearly that due to the “end effects” one cannot build up the poten- 
tial near the end of the electrode, by finite flux elements alone. Rather one has 
to resort to the use of the point source placed at the electrode extremity. Of 
course, even then one does not obtain a strictly uniform potential over the 
electrode surface; however, a potential distribution whose maximum differ- 
ences do not exceed 2 percent of the average potential can be obtained with- 
out very much difficulty and for practical purposes this is of course fully 
sufficient.* 

As a practical aid, tables were constructed giving log I'(y), log sin wy 
and 7 ctn ry for OS yS1 the values of y increasing by intervals of 0.005. 
From these, tables of ®o(x, w) and ®,(x, w) were computed for various values 
of x and these served to give the PA. 











a 





* Although the above analysis has said nothing of the bottom of the electrode it really 
corresponds to an almost hemispherical cap. From Eq. (29) it follows that ®o(x, p2, x) =®o(x, 
0, x+p2)+log 2. The bottom of the electrode is therefore not exactly spherical, but is slightly 
flattened. 
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vi 
ui 


In Fig. 7 are plotted the resultant values of ® over the electrode surface 
and on its geometrical extension to the bottom of the disk, i.e., to w=0.5 for 
various values of x. Their analytical representations are given by the follow- 
ing equations in which are included the complete penetration and zero pene- 
tration cases. The functions ®) and ®, in the equations below are really the 


a 


(Xx, P, W) 
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Fig. 7. Curves giving the potential distributions over the electrode surfaces of radius p2= 
0.001, after adjustment of the flux distribution according to equations shown below. Dotted 
curves represent the unadjusted distributions. 


functions of Eq. (19) and (17a) which reduce to the corresponding functions 

of Eq. (36) and (36b) when p is equal to po, the effective electrode radius. The 

coefficients have been adjusted so that the resultant ® reduces to the uni- 

form value ®, when p2=0.001. The constant ®, has then been added to the 

expressions so as to make ® vanish on the electrode surface. The results are: 
x =0.50: electrode penetration = 100%. 


2 
$(0.50, p, w) = 15.2018 — 2 log — 
p 


x = 0.45: electrode penetration = 90%. 


(0.45, p, w) = 15.62 — {1.4336 (0.45, p, w) — 0.2334 (0.44, p, w) 
— 0.054,(0.42, p, w) — 0.054,(0.40, p, w) — 0.02,(0.375, p, w) 
— 0.044,(0.30, p, w) + 0.00386,(0.45, p, w)}. 


x =0.375: electrode penetration = 75% 


o™ 


(0.375, p, w) = 14.41 — [1.354,(0.375, p, w) — 0.184)(0. 36, p, w) 
— 0.07,(0.34, p, w) — 0.054,(0.30, p, w) — 0.054 (0.20, p, w) 
+ 0.00426,(0.375, p, w) + 0.00026,(0.37, p, w)}. 
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x =0.25: electrode penetration =50%. 


(0.25, p, w) = 13.28 — }1.4540(.25, p, w) — 0.244,(0.24, p, w) 
— 0.104 (0.22, p, w) — 0.0314 (0.20, p, w) — 0.07 40(.15, p, w) 
+ 0.0038®,(0.25, p, w)}. 

x =0.125: electrode penetration = 25%. 

#(.125, p, w) = 11.25 — {1.3440(0.125, p, w) — 0.224 9(0.115, p, w) 
— 0.074,(0.10, p, w) — 0.054 (0.08, p, w) — 0.034 9(0.05, p, w) 
+ 0.00354,(0.125, p, w)}. 

x =0.05: electrode penetration = 10%. 

(0.05, p, w) = 9.68 — 11.254 (0.05, p, w) — 0.194 ,(0.04, p, w) 
— 0.064,(0.02, p, w) + 0.00333,( .05, p, w)}. 

x =0.025: electrode penetration =5%. 

(0.025, p, w) = 8.30 — 31.304 9(.025, p, w) — 0.264,(0.02, p, w) 
— 0.044,(0.01, p, w) + 0.00284,(0.025, p, w)}. 

x = 0.001 : electrode penetration = 0.2%—‘“nonpenetrating electrode”. 

(0.001, p, w) = 2002.308 — ,(0, p, w). 


To illustrate more clearly the effect of the superposition, the original un- 
adjusted potential functions Pp(x, p, w) have been plotted in Fig. 7 for a few 
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Fig. 8. Flux distributions over electrode surfaces required to give the potential distributions of 
Fig. 7. @p®p represents a point source of strength ap. 


values of x for p=p2=0.001. These are the broken curves. From these it ap- 
pears that near the bottom of the disk the potential due to a uniform flux 
along the electrode and that due to a uniform potential over the electrode 
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surface run closely parallel to each other. However, near the electrode ex- 
tremities the potentials fall off much more steeply for the adjusted potential 
distribution. This of course corresponds to the excess flux density near the 
bottom of the electrode which was added to make the potential uniform over 
the electrode surface. 

Since the physical meaning of the superposition method used above is that 
of the superposition of flux elements, a plot of the flux distributions may be 
of interest. These distributions have been plotted in Fig. 8 for p2=0.001, and 
various values of x. As in no case does the flux density exceed 2, the scale has 
been shifted for the various values of x. The discontinuous curves of Fig. 8 
are clearly the approximations implied by the superposition method, of the 
continuous flux distributions which would give a perfectly uniform potential 
over the electrode surface. By drawing smooth curves through the discon- 
tinuous segments, one should get a still closer approximation to the actual 
continuous flux distribution over the electrode surface. 

These curves show that not only does the flux always require a concen- 
trated point sink at the electrode extremity, which must be added to the 
more continuous distribution along the electrode proper, but also this latter 
distribution appears to be too concentrated near the electrode extremity to be 
represented by any very simple continuous function. One may note, too, that, 
as is to be expected, the relative steepness of the rise of flux near the elec- 
trode extremity decreases with increasing penetration. This, of course, cor- 
responds to the decreasing relative importance of the nonradial part of the 
flux with increasing penetration. 

Before closing this discussion it may be of interest to show more clearly 
that the above superposition method is essentially nothing more than the 
direct approximate solution of the integral Eq. (30). Thus, it will be clear 
that the superposition of the functions ®, A®, and ®, simply amounts to 
writing the resultant potential functions, except for a constant, in the form: 


p= Yam Po( Xm, p, Ww) 
which according to earlier formulae may be expanded as: 


p ~ Ko(2nzp) ' 5 
7 = % ———— cos 2nrw} y= sin 2nt7Xm + 2nmdy COS 2nrx py} 
1 nT 


2 
+ (ap + > xXndm) log —, ; 
p 


where a, represents the strength of the point source at x =X». 

Returning now to the integral Eq. (30) it will be noted that the most 
direct analytical approximation method that could be used to solve it would 
be the following. First one would replace the integral fo7a(y) sin 2n7ydy by 
> mdm Sin 2NTXm, Xm <x, with evenly spaced x,, and then after choosing pz 
one would set the right side of Eq. (30) equal to a constant for each of a set 
of values of w, not exceeding in number the number of values of x». These 
would give a linear system of equations for the coefficients a,, which could be 
solved by the usual algebraic processes. 
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Now as seen by comparing Eqs. (30) and (31), a(u)~q’(u). Hence as 
q(u) varies but slightly for small values of w, the a,, for small values of x,, 
would be small and can be neglected to begin with. Further, it appears that 
due to the very rapid drop in ® for w=x anything but a very close spacing 
of the x,, would give considerable oscillations in ® between the x,,. From these 
considerations it is seen that the natural way of treating the above proposed 
linear equations in a, would be to deliberately drop the terms a,, sin 27x, 
for small x,, and then avoid the close spacing near x,,=x by adding a point 
source at x, =x. But this procedure is clearly essentially identical with the 
method actually used and expressed by the above formula for ®. It is seen 
therefore that the present superposition method is really a direct method of 
solving the integral Eq. (30) in which, however, the resulting algebraic equa- 
tions were simplified and solved graphically by applying the available a 
priori information as to the nature of the functions Po(x,,, p, w). 


D. Potential distribution for large values of p with uniform electrode poten- 
tial 
To see how ® behaves when p is large, it should be recalled from the last 
section that for practical purposes ® can be represented by the expansion: 


(x, p, w) = ® — YoanPo(xm, p, W)i xm S 4, (37) 


where ®, is the constant potential over the electrode surface, the Bp (xm, p, w) 
are the potential functions of Eq. (19) for x =x,,—representing the potential 
due to a uniform unit flux density extending from »=0 to p=x,,—and a, 
is an effective “strength” associated with ®o(x,,, p, w) i.e., it gives the numeri- 
cal flux density for the particular flux element extending from nu =0 to p=Xp. 
It is to be understood, further, that at least one of the terms ®o(x,,, p, w) is 
to be the potential function for a point source, usually situated at x,, =x, the 
potential of which is given by Eq. (17a) with gdu=1, u=x. 

Now as already pointed out the analytical equivalents of the functions 
Po(xm, p, W), which are particularly suitable for use at large values of p, are 
given by Eq. (21): this means that: 


; £ 1 2 
Do(Xm, p, W) = 4 |— don —Ko(2nmp) cos 2nrw sin 2urx», + X» log =| . (37a) 
us l n p 


Since, as has been seen in the earlier sections, the terms in the summation 
are very small for values of p21, one has immediately as the approximation 
to By(xm, p, w), valid for large p: 
? 
Po(tm, p, W) = 4x» log —- (38) 
p 
Similarly: 


2 
P, (4m, p, w) = 4 log — - (38a) 
p 
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One can therefore rewrite Eq. (37) in the form: 


) 
m@ = h. — 4 > mtn log - aT 9, 
p 
where x,, = 1 for the a,, referring to ®,. 

It is thus seen that for large values of p, the potential is for all practical! 
purposes independent of the depth and varies logarithmically with the radia! 
distance from the electrode. This is more clearly seen from Fig. 9 in which 
are plotted a number of equipotential curves for x =0.25, p,=0.001, p; =2 
The potential functions used are really the a,,®) of Eq. (36) since the term 
®. adds nothing to the nature of the distribution. 
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Fig. 9. Equipotential curves in a disk with a 50 percent penetrating electrode. is the function 
given by Eq. (37). R=fraction of total potential drop across the disk. 





Again from Eq. (39) one may get a convenient check with regard to the 
total flux passing through the system: As in Eq. (22), this is given by: 


9 


1/2 AO 
Q=- dna f p — dw = 41a ) amXm. (40) 
0 dp 


This result is quite rigorous, since even if the exact Eqs. (37a) had been used 
the trigonometric terms would have vanished on integration. 
E. The resistance of the disks as a function of the penetration 


Recalling that Eq. (39) gives the potential drop between the electrode 
and external boundary, at pi, the resistance of the system can be obtained as 
in Part I by dividing A®, of ® of Eq. (39), by the total flux Q, as given in 


Eq. (40). Thus: 
o.*% 1 [ a | (41) 
Q = — = — | ———_ — log — ]. 
QO Ta 45 mdm . Pi 


For purposes of plotting, however, it is more convenient to use the di- 
mencsionless function 
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~ 


af 1 [ ®, ; =| ae 
= — | ————- —_ log — 41a) 
2 2r 4 XmOm . py 


which gives the resistance of the system for unit thickness of the disk. In 
Eq. (41a) the ®, and 42x,a,, are to be obtained from the expressions in 
Section C, even for the cases x = 0.50 and x =ps. Assuming p1/p2=7;/1r2 = 2000, 
the curves in Fig. 10 were plotted from Eq. (41a). Logarithmic paper was used 
because of the large variation in Qa/2. The difference between the curves I 
and II is clearly a measure of the dependence of the resistance of the disk 
upon its absolute thickness. Thus it appears that the effect of the absolute 
thickness decreases with increasing penetration. 








—— — oo om ——s Basel tA 
400} —-—.—-—- 4 
’ T 
20g-—1— 

















§ 838 








nN w 
Oo 9 
| 
} 
| 
I 





as 
wY WwW Bb OD ®WO 
| | 


a a SA ME ADE 

aNNES MNS mS a ELE 

_ 4 ++} 
| 


totes 


=RESISTANCE OF DISC 





t 
Berle 
} 




















||| / | ii 
OI a2 04 06 102-3 4 6 8102040 60 100 
°*e PENETRATION OF ELECTRODE 

Fig. 10. Resistance of the disk as a function of the electrode penetration. a 2/2 =resistance 

of disk per unit disk thickness. Curve I: p2=0.001, :/p2=2000 with electrode potential uni- 

form, as given by Eq. (41a). Curve II: p2= 1/600; p:/p2= 2000 with electrode potential uniform, 

Curve III: p2=0.001, p:/p2=2000 with uniform flux density over electrode (Eq. (42)); crosses 

represent values computed by Eq. (43). Curve IV: p2=1/500, p:/p:=1000, (computed by Eq. 
(43)). Curve V: p2= 1/300, p:/p2= 1000, (computed by Eq. (43)). 





It will be seen that after a penetration of about 5 percent curves in Fig. 
10 are approximately linear with a slope a little less than one. However, the 
marked curvature appearing for penetrations close to 100 percent is quite 
real and its existence has, in fact, been verified experimentally.* It is of 
interest, too, that the resistances for the “nonpenetrating electrodes” given 
by the 2 points at the extremities of curves I and II fall, without adjust- 


* The writer is indebted to Mr. R. D. Wyckoff and Mr. W. O. Smith for these results 
which were obtained by high frequency experiments performed with a large cylindrical tank 
containing a CuSQ, solution. 
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ment, on the natural extrapolations of the main parts of the curves. There are 
thus no discontinuities or sharp curvatures, so that interpolation for any in- 
termediate points can be made with but little uncertainty. 

Since these computations of the resistance by Eq. (41a) requires the 
preliminary tedious determination of ©, and the constants x, for each elec- 
trode radius or disk thickness, it is of some interest to develop an approximate 
formula that will avoid as much of this labor as possible. Since it is readily 
shown that to abandon the infinite set of images is a far too crude approxi- 
mation, the only other reasonable approximation is to see what can be done 
with a uniform flux density along the electrode, as developed in Part II. 
In this case Eq. (41a) would reduce to: 


a2 1 & , 3 , 
ood bag eo | (41 
, tlhe ; 


However, ®,. which represents the electrode potential must still be de- 
termined. Since the potential over the electrode surface is not constant when 
the flux density is uniform ®, should correspond to an average of the surface 
potential. As a first approximation, however, one may assume ®, to be the 
potential at the top of the electrode, which by Eq. (36) is: 

2r 
®, = $o(x, 0) = 2 log — — 4 log I(x) — 2 log sin rx. 
p2 


Eq. (41b) may then be written as: 


aQ 1 E an 1 (1 — x) , | ” 
— = — | — log — + — log ————_ —- log — ]. (42) 
2 2x L2x . po 2x . r(i + x) . pi 


This is plotted in Fig. 10 for pp =0.001 and p;=2.0 as a function of x, as 
curve III. It will be seen that the aQ/2 of curve III are almost uniformly 
10 percent higher than the correct values given by curve I. Although this 
agreement would be good enough for some purposes it will be presently seen 
that a slight change in the assumption for ®, will improve the approximation 
considerably. But before leaving Eq. (42) it is of interest to compare it with 
the corresponding result of Samsioe who also assumed a uniform flux density 
and took ®, as the potential at the top of the electrode. In the present nota- 
tion Samisoe’s result is equivalent to: 


aQ 1 1 2x 1 1+ 1.44x\°-4 2 
— = — | — log — + — log (=) — log — |]. (42a) 
2 2r L2x p2 2x 1— 1.44% pi 





Since a numerical check easily shows that: 


( + ae fl r(i — x) 

r(i + x) 
to a high degree of accuracy, it is clear that the present treatment gives the 
same result as Samsioe’s for the net resistance of the system under the con- 
dition of uniform flux density along the electrode surface. 
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Returning now to the above assumption that ®.= (x, 0) it is evident 
that it would be more reasonable to choose for ®. the value of ® at some 
definite position down along the electrode. By trial it has been found that 
if ®. is taken as the value of ® at w=3/4x the resistances computed accord- 
ing to Eq. (41b) agrees with the values computed for Eq. (41a) to within less 
than 1/2 percent. The agreement may be seen better from Fig. 10, where the 
crosses are the values computed as stated above, i.e., from the formula: 











12 1 [bo(x, w = 3x, po) 2 
ae ae | : - io aa log =|, (43) 
2 2r 4x Pi 
or, in detail 
aQ -| - les 2x » log F(1.75)F(0.252) 
= —| —< 2 log — — 2 log 75x 25x 
2 2 =e . p2 
‘ (43a) 
— log sin 1.75rx sin 0. 2saxt — log]. | 
y Pi 


By the use of this formula one is able to avoid the tedious work of ad- 
justing flux elements so as to obtain a strictly uniform electrode potential, 
when computing the resistance of the system. Thus with this expression 
curves IV and V of Fig. 10 were computed. These show the variation of aQ/2 
with the electrode radius since they correspond to curves I and II except 
that the electrode radius has been doubled. As the curves show, the effect is 
not large, and as is to be expected aQ/2 varies with the electrode radius 
logarithmically, as for pure radial flow, for large penetrations, and approaches 
an inverse first power variation for the small penetrations. 

It is to be understood that if one is interested in the details of the potential 
distribution the above approximation has no meaning and one must use the 
more lengthy functions of section C. It should be pointed out, also, that since 
the resistance of the system is inversely proportional to the thickness of the 
disk, the resistance for a disk of thickness a/2 will be given by the ordinates 
in Fig. 10, divided by a/2. As to the units it is clear that since aQ is dimension- 
less, the units of 2 will correspond to those of the length a and the specific 
resistance, which has been taken as unity in all the above considerations. 

The above discussion would sufficiently complete the analysis from the 
electrical point of view. However, from the point of view of fluid flow through 
a sand, one is usually interested in the flux or production passing through the 
system for unit potential or pressure difference. That is, the main interest is 
in the reciprocal of the resistance. It will therefore be useful to reinterpret 
some of the above results from this point of view. 

Thus the flux per unit potential difference or production per unit pressure 
difference and per unit sand thickness is given by Eq. (41) as: 
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In Fig. 11 the quantity 2Q/aA® is plotted for several values of ps with 
pi/p: = 2000. For p2— 1/600, 1/1000 Eq. (44) was used. For p2 = 1/400, 1/1600, 
the reciprocal of Eq. (43) was used, i.e.,: 
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Again the differences between the various curves indicate that the flux for 
various penetrations is not strictly proportional to the absolute disk thick- 
ness. Since for a given electrode or well radius the sand thickness is inversely 
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Fig. 11. Curves showing the total flux through the disk for unit potential difference and 
unit disk thickness, for various disk thicknesses as a function of the electrode penetration. p1/p2 
= 2000. Permeability or conductivity =1.0. Curve I: p2=1/400; Curve II: p.=1/600; Curve 
IIT: pp=1/1000; Curve IV: p2=1/1600. Curve V: Flux if the flow were purely radial. 


proportional to pe, the fact that the curves for higher pz lie higher in Fig. 11 
shows that the flux per unit disk thickness is greater for the smaller thick- 
ness, or that the flux for a given percent of penetration does not increase 
quite as fast as the disk thickness. 

The dotted line in Fig. 11 is the flux that would pass through the system 
if the flow into the electrode were completely radial for all penetrations. The 
difference of about 20 percent between this line and the curves above it, 
shows at once that the actual flow system cannot be considered as a super- 
position of a radial flow into the electrode proper and a semispherical flow 
into the electrode extremity, since the contribution due to the latter is only 
about 2.5 percent of the radial flux. Rather, the flux distribution even ap- 
preciably above the electrode extremity is distorted from a radial character 
and many of the flow lines from a depth below that of the electrode come 
up and end on the electrode proper. 

In order that the results of the above analysis be as general as possible, 
no attempt has been made to put in absolute numerical values for any of 
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the dimensional quantities. When such specific results are desired they can 
be readily obtained from those presented by the multiplication or divison 
by the proper coefficient. Thus with regard to Fig. 11, if the actual flux is 
desired for a given pressure differential and disk thickness one need only 
multiply the numerical ordinates of Fig. 11 by the given pressure differential 
and thickness, these being expressed, however, in the same units as the per- 
meability, which again has been taken as unity. 

The writer wishes to thank Mr. R. D. Wyckoff and Dr. P. D. Foote of 
the Gulf Research Laboratory for the opportunity of discussing with them 
questions arising in the preparation of this paper. 
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Chemical Treatment of Rotary Drilling Fluids* 


By H. C. Lawton, H. A. AMBRosE, and A. G. Loomis 
Mellon Institute of Industrial Research and Gulf Research Laboratory, Pittsburgh, Pa. 


(Received March 11, 1932) 


The desirable properties of rotary drilling fluid in the technique of boring oil wells 
are considered. An analysis of the problem of elimination of cuttings by differential 
settling with reclamation of the colloidal content of clay dispersions shows that vis- 
cosity should be lowered as much as possible when the density is held constant. It is 
demonstrated that lowering of viscosity results in a small lowering of pump pressure, 
even though the velocity of the mud-fluid is increased to maintain constant carrying 
capacity for cuttings. Several chemicals have been found that cause remarkable lower- 
ing of the viscosity of clay drilling fluids. Suitable adjustment of the concentration 
of caustic soda in these chemical reagents imparts desirable thixotropic properties to 
the clay dispersion. The theoretical results of this study are applied to some of the 
practical problems of oil-well drilling. 


I. INTRODUCTION 


HERE are two standardized methods now in general use for boring oil 

wells. One method of drilling, known as cable-tool drilling and de- 
scribed as a “standard rig” in oil-field parlance, uses a heavy drilling bit 
attached to a wire cable. By a suitable mechanism, causing withdrawal 
and release of the cable, the bit is alternately raised and dropped, striking 
hammer blows on the rock formations, with the result that the rock beneath 
the bit is crushed and the formation is gradually penetrated. Small amounts 
of water introduced into the bore hole serve to suspend the cuttings and leave 
a more or less fresh rock surface exposed for each blow. The bit is frequently 
withdrawn and a bailer is intermittently introduced to remove the debris 
which, if allowed to accumulate, would soon render the drilling ineffective. 
When a water-bearing formation is encountered and the hole is filled with 
water, the downward velocity of the bit is greatly diminished and drilling be- 
comes relatively slow and inefficient. Such water horizons must, therefore, be 
cased off, entailing considerable expense for pipe and necessarily a reduction 
in size of the bit subsequently used. Thus each water horizon encountered re- 
quires a separate string of casing for sealing-off and the bore hole rapidly 
tapers to an excessively small bore. 

To meet this and other shortcomings of the “standard rig” the “rotary 
drill” is used. It consists of a hollow drill stem carrying a suitable cutting head 
at its extremity. The stem is rotated continuously from the surface, with pres- 
sure applied at the bit. By means of pumps a suitable drilling fluid is circu- 
lated down the hollow drill stem, through orifices at the bit face, and thence 
back to the surface through the annular space between the drill stem and well 
bore. This drilling fluid or mud has three distinct functions, namely: 


* An elaboration of this paper has been submitted by Mr. Lawton as a master’s thesis, 
University of Pittsburgh. 
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(1) It washes the rock face clean and carries cuttings to the surface, where 
the cuttings are settled out in a pit and the mud fluid is used again. 

(2) Water, gas and minor oil-bearing formations encountered during drill - 
ing are sealed off by the hydrostatic head of the fluid column, controlled by 
proper adjustment of the density of the drilling fluid. 

(3) The pore-sealing effect on the face of the bore hole results in the pre- 
vention of caving of loose and unconsolidated formations. Incidentally, the 
colloidal nature of the mud-fluid serves as a lubricant for the drill stem and 
for running casing into the hole. Thus the rotary drilling method allows more 
or less continuous drilling through one or more water or gas horizons without 
necessarily requiring casing of the bore hole. Obviously, if due care is not 
taken, or if improper mud-fluid characteristics obtain, temporary suspension 
of mud-fluid circulation may cause settling of the cuttings, resulting in hope- 
lessly jamming the drill stem and preventing subsequent retrieval. Likewise 
improper control may entail loss of one or more of the important functions of 
the mud-fluid. 

It is apparent that the mud-fluid used is an extremely important item in 
drilling oil wells and it is of interest to mention that the special mud-fluid used 
in certain drilling operations in the past has occasionally represented an ex- 
penditure as high as $30,000; hence, the importance attached to this ap- 
parently trivial item of drilling equipment. It is not too extravagant a stat- 
ment to say that the proper quality and control of the mud-fluid is the key to 
successful drilling. 

It is only recently that petroleum technologists have had a full apprecia- 
tion of the importance of adequate control by chemico-physical means of 
rotary drilling fluid in the technique of boring oil wells. Any drilling fluid may 
be regarded as a system of three main constituents, as follows: water, colloidal 
matter, consisting of both gel-forming and non-gel-forming colloids, and larg- 
er particles suspended in the medium. Of the two types of colloid present, 
the gel-forming type is definitely the more important, for this group imparts 
the necessary properties of an efficient drilling fluid, i. e., stability of suspen- 
sion, viscosity, pore-sealing, and gel-characteristics in general. Colloids of 
the non-gelling type, together with larger particles, contribute somewhat to 
viscosity but more particularly to the density of the drilling fluid. Adequate 
control is a problem of colloid physics and it will be shown how the gel-type 
of colloid present responds to treatment, giving desirable and necessary 
characteristics of drilling fluid, modified, of course, by the inert material al- 
ways present. 

In this paper we shall restrict ourselves to clay drilling fluid and a con- 
sideration of the problem of lowering the viscosity. It is generally desirable in 
drilling practice to maintain the viscosity in the neighborhood of 15 centi- 
poises, as measured by the Stormer viscometer at 10 revolutions per second. 
For a given clay this naturally limits the maximum density of drilling fluid 
that can be mixed without treatment, but the mud may then lack the de- 
sired weight and corrective measures become imperative. A consideration of 
the basic functions required of a drilling fluid, insofar as cuttings and gas are 
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concerned, together with the characteristics that are required to fulfill these 
functions, leads us to the following classification: 
1. Carrying cuttings 
Requires alone or in combination: 

High density | 
High velocity 
High density | 
High viscosity { 


in the turbulent region 


in the viscous region. 


2. Elutriation or dropping of cuttings 3. Release of gas 
Requires alone or in combination: Requires: 
Low density Low viscosity 
Low viscosity Time, or reduction of path 
Low velocity Other factors outside the 
Time, or reduction of path. scope of this paper. 


From the above it is seen that the requirements for carrying of the cut- 
tings up the bore hole and settling of the cuttings in the mud pit are diametri- 
cally opposed; also, one requirement for release of gas is opposed to one of 
the requirements for carrying of cuttings in the viscous region. Low viscosity 
favors two of the desirable functions of the mud-fluid, while a high viscosity 
favors one, the capacity for carrying cuttings in viscous flow. With respect to 
the carrying of cuttings, however, an increase in the mud-stream velocity or 
the density or both will compensate for lowering the viscosity. When other 
consderations, such as ease of handling, reduction of the load on the pumps, 
reclaiming of mud-fluid by elimination of cuttings and gas, etc., are taken into 
account, it is clear that every effort should be made to keep the viscosity of a 
mud-fluid at the lowest possible figure. At the same time the mud weight 
should be maintained at its initial value. The mud could, of course, be diluted 
with water, with resultant lowering of viscosity and elimination of cuttings in 
a mechanical separator; but when the mud is rethickened to the weight re- 
quired, there is danger that the viscosity will be excessive. 


Il. Errect or LOWERED VISCOSITY ON PUMP PRESSURE 


A consideration of the rates of circulation in general drilling practice 
shows that the critical velocity for mud-fluids ordinarily used is exceeded in 
all except the 16 and 20 inch holes and that flow is in the turbulent region in 
the drill pipe in all instances. In the case of turbulent flow Fanning’s formula 
for flow in a circular pipe may be applied to flow in an annular space, such as 
exists when drilling, provided the mean hydraulic radius is used in the for- 
mula.' From Fanning’s formula, p= Kn°?v'-%p°* for unit size of hole of unit 
depth, where # is the pressure drop, K is a constant, 9 is viscosity, v is ve- 
locity and p is density, it is seen that pump pressure will be lowered by only 
the 0.2 power of the drop in viscosity for the same rate of flow and density of 
mud-fluid. As mentioned above, however, the velocity or density or both 
must be increased to maintain the same carrying capacity for cuttings when 


1 Engineering Experiment Station, Univ. of IIl., Bulletin No. 222. 
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the viscosity is decreased. The carrying capacity v, will be given by v,.=v—v,, 
where v is the velocity of flow of the mud-fluid and 2, is the velocity of slip of 
the cuttings. In turbulent flow the velocity of slip between the mud-fluid and 
cutting is given by Rittinger’s formula (for a spherical particle in water): 


d(pi — =) be 





7, = 9 —— - 


: \ p2 
where v, =velocity of slip in feet/sec., d= diameter of round cutting in feet, 
p:=density of cutting, lbs./cu. ft., and p2=density of mud-fluid, Ibs./cu. ft. 
(The constant 9 becomes about 3.5 for flat particles.) 

It is reasonable to assume that this expression is valid so long as the mud- 
fluid flow is in the turbulent region, even though the differential velocity of 
fluid and cutting is sufficient to give a Reynold’s number below the critical 
value. At the maximum the viscosity probably enters the constant (friction 
factor) of Rittinger’s formula to about the 0.2 power or less for liquids show- 
ing anomalous viscosity, such as mud-fluid, approaching a value of zero if the 
drill hole and cuttings are considered to be very rough, as is actually the case. 

From Fanning’s formula, p= K,n°*v!-’ for unit size and depth of hole and 
unit density of mud-fluid. The velocity of flow of the fluid is then given by 
v=(p/K,n°*)'/"-’. The velocity of slip v,=Ke, from Rittinger’s formula, as- 
suming that the flow of fluid past the cutting is turbulent and that the power 
index of the viscosity approaches zero for perfect turbulence. Then 


? = (p/K.n°-?)! 1.8. Ko. 


Solving, p= K.n°?(v.+K_2)!°=K3n°? for constant carrying capacity. 

It is to be concluded, then, that pump pressure varies as n°, approxi- 
mately, for constant carrying capacity of the mud-fluid for cuttings. The 
final result of lowering the viscosity of the mud-fluid will be a small drop in 
pump pressure, while the efficiency of removal of the cuttings remains approxi- 
mately constant. 


Ill. Turxotropy oF Mup-FLvuIpD 


The chief objection to the use of a low viscosity mud in the hole, although 
partially compensated for by increased mud-stream velocity so far as carrying 
capacity is concerned, is the fact that rate of settling of cuttings will be greater 
in case circulation is unavoidably stopped. Rate of settling is independent of 
the height of the column, but the amount settled around the drill stem de- 
pends on the height of the column, the time of settling and the content of 
cuttings. For a given height of column both of the latter factors may be modi- 
fied. A chemical reagent may be used that gels the mud-fluid quickly while 
in a static condition, and settling takes place only during the few seconds 
elapsing before gel-formation: The viscosity of the gel is so large that settling 
ceases and even large particles are held in suspension. When the pumps are 
again started the gel breaks easily with restoration of the original low vis- 
cosity of the mud-fluid. The content of cuttings may be reduced by using as 
high a mud-stream velocity as possible for a given rate of drilling and by 
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settling as many cuttings as possible in the settling pump. It is clear, then, 
that a chemical reagent must lower the viscosity as much as possible, and 
cause the formation of a light gel when the mud-fluid is in a static condition. 
Moreover, the chemical must cause the release of particles larger than 325- 
mesh only and must stabilize by deflocculation particles smaller than 325- 
mesh. In this way muds may be reconditioned and the colloid content re- 
stored. 
IV. GENERAL CONSIDERATION OF THE VISCOSITY FACTOR 

1. Stokes’ law 


From a qualitative consideration of the functions of a drilling mud it has 
been shown that viscosity should be as low as possible for a mud of a given 
weight. The same conclusions follow for dropping cuttings when Stokes’ 
law is considered. 


where V =velocity of fall of a particle, g=acceleration of gravity, r=radius 
of falling particles, p; =density of particles, p2=density of medium, and 7= 
viscosity. In cases where there is no gel formation, hindered settling does not 
begin until 30 percent solids are present in the fluid,? and therefore Stokes’ 
law of unhindered settling applies. For a given weight mud it is clear from the 
formula that the only practical way to increase the rate of removal of cuttings 
is to decrease the viscosity. The same equation with the same conclusion ap- 
plies in the case of release of entrained gas from a mud, since the rising bub- 
bles of gas also follow Stokes’ law approximately. 


2. Differential settling 


In order to reclaim or recondition clay mud-fluids it is necessary to settle 
out the fine sand as well as large cuttings, while at the same time recovering 
the finely-dispersed clay particles. For practical purposes we may say that 
all particles larger than 325-mesh should be settled out while all particles 
smaller than these should be held in suspension. Suppose that we divide the 
mud particles into two main sizes, large ones and small ones, designated by 
Land S, respectively, then we shall have for the velocities of fall of each size 

7 Krz7(pr — pe) “ Krs*(ps — pe) 
L= and Vs = ———————__ 
n n 





We are interested in the practical case that the large particles settle out 
faster, i. e., Vz > Vs, and wish to investigate the magnitude of V; — Vs. 


, 


Vi — Vs = — | r%oe = p2) - rs*(ps _ r») | 


0 


Kri? { m\ re \)\ 
= Efon~ (Z)on~ lt - (=) |} oo 
n 'L rL 


2 Bignell, Oil and Gas Journal, May 22, 1930, p. 50. 
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r,\? r.\? 
n>(2)neo[-(2)] 
'L rL 


This immediately tells us that we should make p, as large as possible and ps, 
p2, and 7 as small as possible. Since p, and ps cannot be subjected to arbitrary 
variation under practical conditions the only conclusion possible is that the 
density and viscosity of the mud should be decreased as much as possible in 
order to obtain a large differential effect. In drilling practice it is important to 
maintain the density as high as possible in order to seal off gas and water 
horizons by the hydrostatic column of the mud-fluid. Our only remaining 
variable, then, is viscosity. If we assume that r,>rs then (75/71)? is negligible 
and 


Vi-—Vs= — (px — p2). 
n 
The final solution again tells us that viscosity must be lowered as much as 
possible in order to reclaim or recondition a mud-fluid by elimination of the 
larger particles by settling. 

The only remaining factor affecting the release of sand and cuttings is the 
length of path of the mud-fluid in the time allowed for release by the capacity 
of the mud storage system. It will be shown that certain chemical reagents 
lower the viscosity markedly; as a result, cuttings settle rapidly and the 
length of the settling ditch ordinarily used is more than sufficient. Baffles 
should be installed in the settling trough in order to break the gel that tends 
to form and the velocity of the mud-fluid stream should be large enough to 
counteract any tendency toward an increase in viscosity through gel-forma- 
tion. The optimum velocity should be just short of turbulence in order to 
effect the most efficient release of cuttings and possible entrained gas. 


V. CHEMICAL REAGENTS 


The equations given above take no account of electrical or adsorption ef- 
fects with which colloid physics deals, for Stokes’ law does not apply to par- 
ticles of colloidal size. By taking advantage of colloid phenomena much better 
results for differential settling are obtained than the above equations would 
indicate. Certain chemicals that are strongly adsorbed by the colloidal clay 
particles cause remarkable lowering in viscosity of the mud-fluid. As a result, 
the large cuttings and sand particles settle out quickly; the small particles 
are stabilized by deflocculation, a light gel is usually formed within a short 
time, and the small particles remain in suspension indefinitely. 

The authors are not aware of any comprehensive quantitative data pub- 
lished on chemical treatment of mud-fluids. Duckham’ refers briefly to the 
use of soda, soda ash and sodium aluminate, with very few data based on 
measurements. Doherty, Gill, and Parsons* mention several chemicals that 
have been used in field practice. A commercial preparation known as “Sta- 


% Duckham, Inst. Pet. Tech., March 1931, p. 153. 
* Doherty, Gill, and Parsons, Oil and Gas Journal, June 11, 1931, p. 52. 
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bilite” appeared on the market about a year ago. Certain private reports 
that have come to our attention recommended sodium silicate as the best 
reagent, but the viscosity was lowered only slightly for the actual mud and 
it was necessary to thin with water in order to eliminate the sand. The same 
results could be obtained without chemical by diluting a little further with 
water. 


VI. LOWERING OF VISCOSITY BY CHEMICAL REAGENTS 


In our experimental work, apparent viscosities of the treated mud-fluids 
were measured with a Stormer viscometer, in a manner previously de- 
scribed.® All samples were screened through a 70-mesh screen before making 
viscosity measurements, for larger particles interfere with the free rotation 
of the drum of the instrument. Practically all the measurements were made 
with Pierce Junction rather than with Seminole mud-fluid, for the initial 
viscosity of the former is considerably higher. The weights of these mud-fluids 
were 9.75#/gal. and 10.31#/gal., respectively, after screening through 70- 
mesh. 

Fig. 1 shows the lowering in viscosity of both Pierce Junction and Semi- 
nole mud when treated with caustic soda. It is seen that concentrations (ex- 
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Fig. 1. Clay mud fluids treated with caustic soda. Viscosity in centipoises vs. percentage caustic 
soda by weight of solids. Temperature 25°C. 


pressed in percentage of caustic by weight of total solids of the mud) of caus- 
tic soda between 0.1 and 0.3 have a marked lowering effect on the viscosity 
of Pierce Junction mud, but that the effect is not so large for Seminole mud. 
Increasing concentrations cause a rapid rise in viscosity and the mud be- 
comes plastic. At still greater concentrations the viscosity drops and the mud 
is flocculated, settling rapidly as a sludge. We would expect the hydroxide 
ion to be the most effective of the monovalent ions because of its strong ad- 
sorption. Further work must be done to compare with multivalent anions. 
Fig. 2 summarizes the data for caustic soda, “Mud-It”, (a commercial pro- 


5 Ambrose and Loomis, Physics 1, 129 (1931). 
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Fig. 3. Pierce Junction mud-fluid treated with various chemical reagents. Percentage change 
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Fig. 2. 


Fig. 3 shows the percentage lowering of viscosity as a function of the con- 
centration of chemical. 
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duct containing sodium aluminate with excess caustic soda), sodium oleate 
(soap), Stabilite, sap-brown (a commercial product containing humates 
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Pierce Junction mud-fluid treated with various chemical reagents. Viscosity vs. 


and lignates), and sap-brown mixed with caustic soda. It is seen that Mud- 
It and sodium oleate have little merit in lowering viscosity while sap-brown 
and Stabilite are superior to caustic soda both with respect to lowering the 
viscosity and elimination of the rapid increase in viscosity shown by caustic 
soda with increase in concentration. 


of viscosity vs. percentage reagent by weight of total solids. 
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Fig. 4 presents comparative results obtained on a different sample of 
Pierce Junction mud with sap-brown, sap-brown plus caustic soda, and Sta- 
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Fig. 4. Pierce Junction mud-fluid treated with various chemical reagents. Viscosity in centi- 
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Fig. 5. Pierce Junction mud-fluid treated with various chemical reagents. Viscosity 
in centipoises vs. percentage reagent by weight of total solids. 
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bilite. The two samples of mud were taken in the field at different times and 
showed differences in color and in initial viscosity. 
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Fig. 5 compares the efficiency of humic acid, humic acid plus varying 
amounts of caustic soda, and Stabilite. 
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Fig. 7. Pierce Junction mud-fluid treated with various chemical reagents. Viscosity in centi- 
poises vs. percentage reagent by weight of total solids. Temperature 25°C. 


Fig. 6 illustrates the action of Stabilite with soda black liquor, both 
alone and with addition of caustic soda. The soda black liquor is a by-product 
of the paper industry, and contains mostly lignins and carbohydrates with 
sodium carbonate. 
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Fig. 8. Pierce Junction mud-fluid treated with various chemical reagents. Viscosity in centi- 
poises vs. percentage reagent by weight of total solids. Temperature 25°C. 


Fig. 7 shows that tannic acid alone has only a slight effect in lowering 
viscosity, while tannic acid with increasing amounts of caustic soda added is 
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superior to Stabilite. Peat with caustic soda also gives excellent results; 
peat contains humic acid chiefly. Sulfite liquor, also a by-product of the paper 
industry, gives excellent results, but somewhat greater concentrations are 
necessary. 

Fig. 8 compares Stabilite with chestnut wood extract, containing tan- 
nates and carbohydrates, and spent liquor; the latter contains the material 
extracted from rice hulls by caustic soda and is a by-product of the Lake 
Charles Products Company. 

From these data it is evident that the organic materials used are much 
superior to the inorganic reagents for the purpose of lowering the viscosity of 
clay muds. The degree of gelling the mud is almost entirely dependent on the 
amount of caustic soda used with a given organic chemical. The mixtures re- 
ferred to in the figures form light gels with the clay. 

It should be emphasized that the water used in these experiments was of 
very low mineral content. Similar results in field practice cannot be expected 
if brines are employed. 
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Radon Content of Soil Gas 
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An electroscope and technique for measuring the radon content of soil gas are de- 
scribed. Measurements were made toanswer several questions regarding the radon con- 
tent of soil gas. It was found that there is no systematic increase of radon content with 
depth of penetration into a soil, that soil formations appear to have more or less char- 
acteristic radon contents, that high radon content of a soil gas is not necessarily an 
indication of the presence of a fault, and vice versa; and finally it is found that there 
is great local variation of radon content of soil gases. 


INTRODUCTION 


N RECENT years there has been some discussion as to the applicability 

of radioactive methods in geophysical prospecting. Various claims have 
been made by different investigators as to the geological significance of the 
measurement of the radon content of soils. It has been stated for instance 
that a fault will be evidenced by a higher radon concentration in the overly- 
ing soil,! and that the radon content of soil in an oil bearing region is higher.” 
In order to test the usefulness of radioactivity measurements for geological 
purposes an extensive study was made of the radon content of soils over re- 
gions whose geology and lithology were rather well known. Measurements 
were made over faulted areas, salt domes, and areas which were geologically 
monotonous. 

The measurements were made in soils which were developed over sedi- 
mentary formations consisting of predominantly sandy clays with some 
slightly consolidated sands, all of tertiary age, by climates with annual rain- 
fall varying from 31 to 50 inches, and average temperature of about 68 de- 
grees Fahrenheit. Most of the measurements were made in wooded areas, but 
some were in cultivated fields and some in natural prairies. 

In a few cases the measurements were from the bed rock below the soil, 
but even in these cases the sampling indicated that the formation had been 
effected by oxidation. With these few exceptions, all measurements are from 
the “B” layer of the soil and are above the level of ground water. 


DESCRIPTION OF APPARATUS AND TECHNIQUE 


The instrument used was a modified form of the Ambronn’ electroscope. 
The technique of sampling and measuring the soil gas which was developed 
after a preliminary study was as follows. 


1 F. Miiller, Zeits. f. Geoph. 3, 330-6 (1927). 
2 L. N. Bogoyavlensky, Bull. Inst. Pract. Geoph. No. 3, pp. 113-124. Leningrad, 1927. 


* Ambronn, Methods of Applied Geophysics. Elements of Geophysics, English translation 
by Cobb, p. 119. 
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A pointed steel driving rod about 1-1/8 inches in diameter was driven 
vertically into the ground to a depth of about two feet. This was then with- 
drawn and an earth auger 1—1/8 inches in diameter was used to continue the 
hole to the desired depth. The purpose of the earth auger was to prevent com- 
pacting the soil and thus reducing the porosity in the immediate vicinity of 
the hole. Two depths of hole were used in the surveys: 3 ft. 9 in. and 5 feet. 
The sampling tube which was then pushed into this hole consisted of a 
length of 5/8 inch iron pipe containing perforations near the lower end to 
allow passage of the soil gas. The lower end of the tube was plugged to pre- 
vent it from becoming stopped up when being pushed into the hole. Slightly 
above the middle of this pipe was welded an 8 inch length of 1—1/4 inch iron 
pipe. This was to provide a seal against the side of the hole when the sampling 
tube was in position. In order to reinforce this seal a spiral of 1/4 inch square 
rod was welded to the outside of this larger section. A hose connection was 
fastened to the top of the sampling tube. The gas sample was thus pumped 
from near the bottom of the hole, through a rubber tube into the ionization 
chamber of the electroscope, by means of a hand operated piston pump. 

The electroscope consisted of a single fiber electrometer mounted on a 
cylindrical ionization chamber of about one liter capacity. The lower end of 
this ionization chamber consisted of a movable piston, thus simplifying the 
emptying and filling of the electroscope. There were two ports in the ioniza- 
tion chamber, one, an exhaust and the other an intake port. In passing from 
the soil to the electroscope the gas passed through two CaCl. drying tubes. 

The procedure in sampling and measuring was as follows: The piston of 
the electroscope was pushed up, thus emptying the ionization chamber. Be- 
fore connecting the hose from the sampling tube to the ionization chamber a 
few strokes of the pump sufficed to replace the air in the system with ground 
air. The hose was then connected to the intake port, the exhaust port closed 
and air pumped from the hole until the ionization chamber was completely 
filled. This method of filling by pumping into the ionization chamber assured 
that the leakage in the system, if any, would be outward, thus preventing the 
dilution of the sample with atmospheric air, which might very easily occur if 
the sample were taken by the suction produced by pulling down on the ioni- 
zation chamber piston. As soon as the electroscope was filled the intake port 
was closed and a stopwatch snapped. Exactly five minutes later the rate of 
movement of the fiber across a given portion of the scale was read. Because 
of the nonlinearity of the scale the same portion was always used in making 
the readings. 

Since in this method, not only the radon but also the active deposit re- 
sulting from the decay of radon, was measured, the radioactivity was rapidly 
changing during the first few minutes after filling the ionization chamber. 
This change (increase) in activity reaches a maximum after about four hours; 
however, it reaches about sixty percent of the maximum in the first five min- 
utes. For this reason the reading was made just five minutes after filling the 
ionization chamber. An additional reason for waiting five minutes before 
reading the fiber movement was the elimination of the effect of the possible 
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presence of thoron. Thoron, being the decay product of thorium X is fre- 
quently present in soil gas. Since thoron has a half life of only 54 seconds, it 
will have decayed to about 3 percent of its initial value in five minutes and 
consequently will have little effect upon the radon measurement. If the read- 
ing is made less than five minutes after filling the ionization chamber and 
thoron is present, its presence and varying concentration during the reading 
will result in an apparently high radon content which will be considerably in 
error. A reading of the blank leak was made immediately before each meas- 
urement in order to be sure that there was no undue contamination of the 
electroscope from the preceding samples. Since a measurement was made ap- 
proximately every fifteen minutes, some contamination was always carried 
over from the preceding readings. This was evidenced by a gradually increas- 
ing blank leakage during the course of the day. The over night rest of the in- 
strument allowed most of the residual active deposit to decay and the blank 
in the morning was always much improved over that of the previous evening. 

The results in this paper are expressed in scale divisions per second of 
movement of the electrometer fiber K 1000. Calibration showed this to cor- 
respond to a radon content of the gas of 15.39 X10~-" curies of radon. To re- 
duce this to a value for the radium content of the soil is difficult because the 
quantity of radium is very small and some of it is probably inside the grains of 
the soil, hence it is not certain that all of the radon escapes to the soil gas 
(according to Satterly, only one sixth to one twentieth of the radon formed in 
the soil gets into the soil gas). Also differences in porosity of the soil will 
cause different holes to draw gas from greater or less volumes of soil. In some 
cases the porosity was found so low that no gas could be extracted with the 
pump. It is therefore only an approximation which can be made, but it is con- 
cluded that 100 on our electrometer scale in this paper will represent about 
0.57 X10-" curies of radon per cc of soil. This electrometer scale which is here 
called “electroscope reading” is the rate of movement of the electrometer 


TABLE I.4 Comparison of values for radon content of soil gases. 

















a . curies 
Author Soil type Content in = 
cm? 
A. Gockel : 0.07 to 0.28 K10-” 
J. Satterly 0.25 xX10-” 
J. Satterly Swamp 0.310-" 
. °. . 94) 

J. C. Anderson Sandstone oO 10-8 g Th 
K. Kahler Sand 0.008 x 10-" curies. 
J. Olujie Loam 1.20 X10-" curies. 
P. Ludewig and E. Uranium containing rock 0.3 to 18.2 X10-# 





Lorenzer 





fiber in divisions per second, multiplied by 1000. Consequently an electro- 
scope reading of 100 corresponds to 1539X10-" curies of radon in the ioni- 
zation chamber. The above figure of 0.57 X10-" curies per cc of soil is ob- 
tained by assuming a porosity of 33 percent. If we assume the density of the 


* “Radioaktivitaét”, St. Meyer and Schweidler, p. 593, 
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soil grains as 2.6 we obtain 0.33 X10-" curies/gm of mineral for an electro- 
scope reading of 100. If we assume that about one sixth of the radon formed 
gets into the soil gas we obtain a value of 1.98 X10-" gm radium per gram of 
mineral. The average value for the radium content of sedimentary rocks is 
1.4X10-" gm/gm. 

Fletcher found by analysis, a radium content in two soil samples from 
Dublin, of 1.3 and 2.9X10-" gm/gm.> Table I shows for comparison the 
values of radon content of soil gases as obtained by different authors. 


RESULTS 


Data on the following points were obtained during the survey: (1) The 
depth at which the radon content is representative of the formation. (2) Com- 
parison of results from areas of different soil texture. (3) Effect of variation in 
moisture content on readings. (4) Whether or not a given formation has a 
characteristic radon content. (5) Whether or not an outcropping fault is char- 
acterized by a higher radon content. (6) Whether the radon content is cor- 
related with any other property of the soil. 


1. The depth at which the radon content is characteristic of the formation 


A survey was made over an area in which the soil was fairly uniform in 
texture: clay at the surface with some sand toward the bottom. Two holes 
were made at each station, a shallow hole, 3 ft. 9 in. and a deep hole, 5 ft. 
deep. The deep hole was placed about one foot from the shallow hole. The 
results given in Table II show that for the depths at which samples were 
taken, there is no correlation with depth. 


TABLE II. Radon measurements from deep and shallow holes. Electroscope reading (division per 
second X 1000). 


Shallow Deep Shallow Deep 
1 0 102 89 
121 131 0 156 
78 66 168 123 
86 58 31 2 
158 131 0 0 
0 0 


It has been stated by other investigators’ that there is an increase in radon 
content with depth of penetration into a soil. In none of our surveys did we 


TABLE III. Radon measurements from deep and shallow holes. Electroscope reading (divisions per 
second X 1000). 


Shallow Deep Shallow Deep 
142 150 49 142 
132 0 252 170 
101 57 241 215 
108 252 136 78 
134 166 106 89 





76 87 61 151 








5 A. Gockel, Die Radioaktivitait von Boden and Quellen, p. 36. 
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find evidence confirming this statement. Table III shows the data from an- 
other survey in an entirly different region than that shown in Table II. 

Again we notice that there is no correlation with depth. 

It was noticed that in regions where the soil near the surface was sandy 
the readings from shallow holes were in general lower than those from deep 
holes. This was due to one or both of two causes: the sandy nature of the soil 
may have allowed the gas to be diluted by atmospheric air through soil 
breathing or the seal at the sampling tube may not have been sufficiently 
tight to prevent leakage. 

In order to determine whether there really was leakage past the seal a 
method was developed for testing each hole. This consisted in introducing 
acetylene gas around the shoulder of the sampling tube just below the surface 
of the ground. The acetylene was prepared from a little calcium carbide and 
water. A U.S. Bureau of Mines carbon monoxide detector (without the char- 
coal adsorbent) was applied at the outlet end of the pump. If this detector, 
capable of detecting a very small fraction of a percent of C2He, showed any 
discoloration the seal was obviously leaking and the hole was abandoned. 
This very sensitive method of detecting leaks showed that a large percentage 
of shallow holes in the more sandy type of soil, leaked and consequently gave 
a lower radon content than the true value. It would thus appear that where 
the soil is dense or where there is a dense layer at or near the surface the radon 
content near the surface will be characteristic of the place from which the 
sample is taken. In more sandy regions the higher porosity of the soil may per- 
mit dilution of the sample with atmospheric air. It will be shown later that 
there is great local variation of radon content and where the deeper gave 
lower readings than the shallow holes this is probably due to such local varia- 
tion. In the Reynosa gravel in those places where it was possible to drill a 
hole the readings were zero. It is impossible to say whether this was due to 
the composition of the gravel or to leakage through the large pores since at 
that time the method of testing for leaks had not been developed. 


2. Effect of soil texture on radon content of soil gas 


Measurements were made on soils ranging in texture from coarse gravel 
to sticky gumbo clay. In coarse sandy soil both high and low readings were 
obtained, and the same was true for the clays. In those regions where a sandy 
soil was encountered below a denser more clayey layer the radon content was 
always quite high. This is probably due to the sealing action of the denser 
clay layer which prevents diffusion to the atmosphere.* However, in general, 
it may be stated that the texture of the soil has little if anything to do with 
the radon content of the soil gas. 


3. Effect of variation of moisture on radon content 


A survey was made in January after a protracted rain over an area which 
was quite low and swampy. In March of the same year after much drier 
weather the same area was surveyed again. In some cases the sampling tube 


6 E. Blanck, Handbuch der Bodenlehre (Springer, 1930), p. 304, Vol. VI. 
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was put into a hole of the original survey. The comparison of results from 
these two surveys is shown in Table IV. 


TABLE IV. Comparison of radon measurements in dry and wet weather. 











Electroscope reading Electroscope reading | Electroscope reading Electroscope reading 














First survey Second survey First survey Second survey 
(div. /sec. X 1000) (div./sec. X 1000) (div. /sec. X 1000) (div. /sec. X 1000) 

51 15 406 7 
64 115 195 23 

115 134 195 45* 
98 93 195 $¢7" 

177 85 208 80 

177 88 28 $7" 

406 0 








* These readings taken by putting the sampling tube down the original holes made on the 
first survey. 

It will be noted that with but three exceptions the results of the second 
survey are all lower than those of the first. This tends to confirm the hypo- 
thesis that in times of heavy rainfall the radon is apt to be sealed in the 
ground and consequently to come more nearly to equilibrium with the con- 
tained radium salts. 


4. Does a geological formation have a characteristic radon content? 


A number of profiles was run across several formation contacts. It was 
noticed in several cases that as the formation contact was approached from 
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Fig. 1. Radon survey over four soil horizons. Numbers indicate electroscope readings (div. 
per second < 1000). Scale 1:75,000. 


the upper part of the formation the radon content increased. The formations 
over which the survey was made are the Wilcox, Carrizo, Mt. Selman and 
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Cook Mountain. The formation contacts were those between the formations 
in the order named above. The accompanying map, Fig. 1, shows the radon 
readings across these formations and contacts. 

The averages of the readings in each formation are as follows: 


Wilcox 34 Mt. Selman 56 
Carrizo 30 Cook Mountain 82 


It would seem from this that although individual readings varied between 
very wide limits the average over a formation may have some significance. 
The dip of the beds was to the east southeast; the lower part of the surface 
of each bed is toward the west in Fig. 1. It will be seen that the radon content 
of the soil near the bottom of the bed seems in general to be higher than that 
in the upper part. This is especially noticeable on another survey made across 
the contact between the Wilcox and the Midway (Fig. 2). 
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Fig. 2. Radon survey over three soil horizons. Numbers indicate electroscope readings (div. 
per second X 1000). Scale 1:63,100. 


There was less variation in the Cook Mountain than in any other forma- 
tion. The Carrizo, a non-marine formation showed the lowest radon content 
of any horizon studied. It is composed of coarse sand, generally very loose. It 
contains cross bedded sandstones at intervals and high readings were found 
near and in these sandstones. The Wilcox in general is quite variable as far 
as radon content is concerned. Rapid changes have been noted by the geolo- 
gist in the texture of this formation, and tests of surface water and surface 
specific resistivity have confirmed this variability. 

In general the results of the survey indicate that formations may have 
characteristic radon contents. Additional surveys over other formations 
would be desirable to obtain further data on this point. 
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5. Does an increase in radon content indicate the presence of a fault? 


It is stated by Ambronn and others that certain “Graben” and faults in 
Europe are detectable by the large increase in radon content of the gas in the 
soil above them. In order to test this statement, fairly extensive surveys were 
made over areas whose subsurface geology was quite well known. Several 
parallel profiles were run across a fault whose surface trace was known to 
within a few feet. These are shown in Fig. 3. Although three of the profiles do 
show a high reading at the true fault line, A, the data would also indicate a 
fault at B. However, it is known that there is no fault at B; consequently it 
is seen that the data are very misleading. 


























Fig. 3. Radon surveys across a known fault. Vertical scale, electroscope readings. Hori- 


zontal scale 1’’=1400’. A—A, position of known fault. B, erroneous fault which would be 
indicated by radon survey. 


Additional measurements were made at the outcrop of a very large fault. 
Measurements made in the soil twenty-five and fifty feet from this fault gave 
readings of 30 and 122 respectively. A sample of gas was taken from the very 
fault plane itself and gave a reading of 99. It is evident that there is no en- 
richment in radon resulting from circulation up the fault plane in this case. 
In regions of Europe where similar measurements have been made the out- 
cropping faults were indicated by higher radon concentrations. This is due 
no doubt to the fact that upward circulation of water along those fault planes 
has carried radon and other radioactive substances to the surface and pro- 
duced the observed effects. It seems safe to say, however, that unless there is 
upward circulation of water to the surface through a fault plane, there will 
be no increased radon content in the superposed soil. 
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In order to study further the extent to which the radon content of a soil 
is the result of extremely local conditions a survey consisting of 16 deep and 
12 shallow holes was made over an area’30 X35 feet.’ As seen’in Fig. 4, the 
holes were arranged in the form of two grids, one superposed on the other 
and so placed that the holes of one grid were equidistant from the holes of 
the other. The soil here seemed to be quite uniform; the surface was entirely 
covered by sod. In drilling the holes a light sandy clay of a bright reddish 
color was first encountered. The soil at the bottom of the holes was lighter in 
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Fig. 4. Radon content of soil over limited area. Figures indicate electroscope readings. XK = 
holes 5 feet deep O =holes 3 feet 9 inches deep. Outside dimensions of rectangle 30’ x35’. 


color, varying between a reddish and a yellowish orange. The bottoms of the 
holes were very sandy. 

As will be seen from the figure the variation in radon content among these 
holes was as great as was encountered anywhere else. In other words, ordi- 
nary local variations in soil will be sufficient to mask any other variations of 
greater geological significance. In attempting to determine a characteristic 
radon value for a given formation a great many measurements must be made 
in order to average out any local effects. 


6. Is radon content correlated with any other property of the soil? 


The measurement of the radon content of the soil is really the measure- 
ment of something quite different. It is a measurement indirectly of the uran- 
ium content of the soil. Radon is a decay product of radium which in turn is 
a decay product of uranium. Radon has only a very brief life, any quantity 
being practically completely decayed in forty days. Consequently it must 
be constantly renewed in the soil gas. It is renewed by the decay of radium 
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which is a solid element and of much longer life, since a given quantity of 
radium will be almost entirely decayed in 16,000 years. The soil at any ap- 
preciable depth has probably been more or less undisturbed for much longer 
periods than that so the ultimate source of the radon is the primary parent 
element uranium which has a life of billions of years. 

The uranium and radium being solids are contained in the grains of the 
soil and consequently only a comparatively small fraction of the total radon 
produced is able to escape from the grains and get into the soil air. Hence the 
radon content of a soil gas will depend not only on the parent radioactive ma- 
terial contained in the soil grains, but also upon the size and nature of the 
surface of the grains. Since the radon is indirectly a measure of the uranium 
present it appears that there should be some correlation, other conditions re- 
maining constant, between the radon and certain heavy mineral constituents 
of the soil, such as zircons and biotites. 


CONCLUSIONS 

A study of the radon content of soil gases over extended areas indicates 
that there is no correlation between radon content and depth, that the radon 
content of a given formation may have a characteristic value, although con- 
siderably more data on this point would be desirable. It is also indicated that 
abnormally high radon content of soil gas does not necessarily show the 
presence of a fault, and that there are very great local variations in radon 
content of such a nature as to suggest that there may be some correlation be- 
tween radon content and certain heavy mineral constituents of the soil. 

The authors wish to acknowledge the help and cooperation of Dr. Paul 


D. Foote, Dr. A. E. Ruark, and Mr. R. L. Chenault in carrying out this in- 
vestigation. 
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Speed, Speed Factor and Power Input of Different Designs of 
Diffusion Pumps, and Remarks on Measurements of Speed 


By T. L. Ho* 
Norman Bridge Laboratory of Physics, California Institute of Technology 


(Received March 5, 1932) 


The speed factor of a diffusion pump is defined as the ratio of its real speed to the 
ideal maximum speed. The fact that the speed factor is much less than unity is due to 
three causes, namely, the vapor dispersed through the slit, the nozzle chamber resis- 
tance and the diffuse reflection of gas molecules at the boundary of the vapor-stream 
and at the condensing wall. Some effectson the speed are considered as being independ- 
ent of the configuration. The following general conclusions result from the measure- 
ment of the speeds of different configurations. (a) In order to reduce the amount of dis- 
persed vapor by the slit, the vapor nozzle should be divergent instead of being parallel; 
(b) in order to reduce the chamber resistance the space surrounding the slit should be 
large and the nozzle chamber should be cold; (c) there exists for each design a certain 
slit area at which the speed is maximum, but there has been no general rule to prede- 
termine it ; (d) in designing a large pump either a narrow throat annular nozzle or mul- 
tiple nozzles can be used. In measuring speeds some precautions should not be over- 
looked, unless the results will be falsified: Either the determined speed will be too small 
or too large (pseudospeed). 


I. INTRODUCTION 


LL present constructions of diffusion pumps can be classified as one 
single or a combination of the three following general types. A, B or C, 
in Fig. 1, show a longitudinal section of the principal part of each type 
which characterizes its inverted and titled forms. One cannot say generally 
that one type is more efficient than the other and that pumps of the same 
type are equally efficient. Each type can have different designs, hence the 
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Fig. 1. The principal parts of three types of diffusion pump. The arrows denote the direc- 
tions of the vapor, LV the low vacuum, HV the high vacuum and j the cooling system. 
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efficiency is different for different individual pumps. Very few previous in- 
vestigators tried to explain the reasons of the efficiency of their designs. How- 
ever, the need of the high speed pumps in modern research renders a discus- 
sion of the physical principles involved in the design of efficient pumps a 
necessity. The present investigation is to determine experimentally the speed, 
speed factor and power input of different designs and to draw the general con- 
clusions by discussing the experimental results. 


* Research Fellow of the China Foundation for the Promotion of Education and Culture. 
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II. DEFINITION OF SPEED FACTOR 


The volume of a gas escaping per unit time, through a slit of area A on 
a thin wall of a gas container, into a perfect vacuum or into a space in which 
the escaping gas molecules are completely and continually washed away by 
some arrangement, is approximately: 


RT \}!2 
y -( ) a (1) 
24M 


where R is the gas constant, 7 the absolute temperature of the gas and M 
the molecular weight of the gas, obtained by considering the number of gas 
molecules striking a surface A per unit time. This is also the ideal maximum 
speed of a diffusion pump having a slit of area A. For air at room tempera- 
ture, V =11.7 liters per cm’. We define the speed factor as the ratio of the 
actual speed of a diffusion pump to the ideal maximum speed. This speed 
factor is of the utmost importance in pump design. 


III. GENERAL CONSIDERATIONS 


Considering a gas passing through a slit into a vapor stream as shown in 
Fig. 2, in the space above the slit, (corresponding to the nozzle chamber) 
according to its different states of motion the vapor is to be considered as 
consisting of two types: directional vapor and nondirectional vapor. The 
directional vapor is directly radiated in from the slit and has a resultant up- 


high 
vacuum 
=— 9% 


odon 


2] °° le 
“7 
low vac. 


Fig. 2. The principal part of a diffusion pump of type A inverted. 


ward momentum, while the nondirectional vapor is produced by re-evapora- 
tion and collision and possesses thermal agitation which depends on the tem- 
perature of the chamber. We first consider that the gas molecules diffuse 
through the directional vapor alone. Gaede' and more recently Molthan? 
treated this kind of diffusion process analytically. They reached the con- 
clusion that the number of gas molecules striking the slit area A per unit 
time is 

z = jaNQA (2) 


where N and Q are respectively the concentration and the mean velocity 
of the gas molecules; and a=f(\/d) <1, where \ = mean free path of gas mole- 
cules in vapor and d=width of slit, which increases with A/d. However, N 
cannot be the same throughout the space above the slit whenever there is a 
flow of gas. So, we regard this N only as the molecular concentration of the 
gas near the slit. Different molecular concentrations along the path of the 


1 Gaede, Ann. d. Physik. 46, 357 (1915). 
2 Molthan, Zeits. f. Tech. Physik. 7, 377 (1926). 
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gas in the chamber can be considered as due to the resultant resistance 
offered by the chamber and the vapor in it, mostly the nondirectional vapor. 
Before discussing this resistance we have to derive another coefficient. 

If there were no chances to return for gas molecules entering the slit, 
there would be z gas molecules carried away per unit time by the vapor 
stream. Actually, there are some gas molecules which rebound from the 
vapor molecules and also from the condensing wall below the slit. 

Let » be the number of the gas molecules returning through the slit per 
unit time, and assume this is directly proportional to z, i.e., 2 =kz. Then the 
total number of the gas molecules extracted per unit time is 


Z=2-N 


4 
tayANQ, 4) 


II 


where y =1—k, <1, and depends largely on the shape of the nozzle and the 
intensity distribution of the vapor stream below the slit. 

To express the speed of the pump, substituting in (4) N=p/m and Q 
=(8p/mp)'?, where m, p, p are the mass of a molecule, the pressure, and the 
density of the gas respectively, and using the perfect gas law, we obtain the 
volume of the gas extracted per unit time, 


§ ( —\- (5) 
S = ayA\| —— , 
2nM 


This is the speed at a point near the slit. 

Now we come again to the chamber resistance. There is a speed difference 
between the inlet mouth and the region near the slit due to the chamber 
resistance. 

Let S, be the real speed (i.e., the speed at the inlet of the pump) and R, 
the chamber resistance. Then, by the definition 3 


t 4 ‘in 
, ” 
or 
; S wi 
"1+ SR, 


But R.=1/U, U being the conductance of the chamber. Regarding a 
rarefied gas passing through a vapor as a gas permeating a porous substance, 
U can be written as: 

U=(RT/M)*'*KF (dimension of chamber, density of vapor), 
where K is a numerical constant, and F is increasing with increasing cross- 
section of the chamber and with decreasing vapor density. 

Then we have, by substituting the values of S and R, in (6), 


. , ( Ay . 

S, = aByA\| —— ] 
2nM/) ’ 

3’ Dunoyer, La Technique Du Vide, p. 132 (1924). 
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where 


1 + ayA ————_ 
2r)'!*KF 
which increases with F. 

By comparing (1) and (7), the product of the coefficients aBy is the 
speed factor of a pump of given area A. Hence in designing a pump of large 
speed one should make each of the coefficients a, 8, y as large as possible. 

In the above expression of speed, the assumption is made that the in- 
tensity of the vapor stream is so large and the fore pressure is so low that 
there are only very few gas molecules diffusing back from the fore vacuum 
to the slit. 

IV. GENERAL EXPERIMENTAL RESULTS 


Some effects on the speed, which are independent of the design of the 
pump, have been experimentally studied by the author, of which only a 
brief summary shall be given. . 

(a) The speed is independent of the fore pressure when the latter is low, 
whereas the upper limit of the fore pressure is quite critical. The same result 
was shown by Langmuir‘ and can be seen also in the speed-pressure graphs 
obtained by Hickman and Sanford.® This critical fore pressure increases with 
the vapor pressure in the boiler. 

(b) If the fore pressure is kept below its critical value, there exists an 
optimum value for the vapor pressure, above which the speed decreases 
again, as stated first by Gaede. This is due to the fact that a and 8 increase 
with decreasing vapor pressure and the weak vapor stream causes an in- 
crease of gas molecules reflected on it and a back diffusion from the fore 
vacuum side. 

(c) At a given optimum range of vapor pressure in the boiler the speed 
is nearly inversely proportional to the square root of the molecular weight 
of the gas to be extracted. This shows, from (7), that the speed factor is 
nearly the same for different gases. Indeed, for a given width of the slit, 
Gaede’s expression of @ is a function of the mean free path, A, of the gas - 
molecules in the vapor, which depends largely on the diameter and velocity 
of the gas molecules. However, when J is larger than the width of the slit, 
a is almost independent of \ and it is the same for different gases. 

(d) The speed is independent of the pressure of the gas being extracted. 
This is to be expected, since (7) does not involve p. Of course, the range of 
the pressure in which the speed is constant is limited by the boundary con- 
ditions of (7). Gaede also experimentally verified it by his observations on 
the first diffusion pump. 

(e) Varying the temperature of the cooling system, the speed first in- 
creases with decreasing temperature and then reaches constancy, since the 


4 Langmuir, Gen. Elec. Rev. 19, 1060 (1916). 
5 Hickman and Sanford, R.S.I. 1, 140 (1930). 
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coefficient 8 is obviously increased by lowering the temperature of the cham- 
ber, and also the rapid condensation of the vapor will increase the coefficients 
a and y. 

(f) The speed of the pump filled with n-butyl phthalate is higher than that 
of the same pump filled with mercury. For the n-butyl phthalate has a lower 
vapor pressure than mercury as has been determined by Hickman and San- 
ford,® hence it causes an increase of a and 8. 


V. EXPERIMENTAL RESULTS OF DIFFERENT DESIGNS 


The word “design” applies only to the size and shape of the principal part 
(i.e., the vapor nozzle and its surroundings) of the pump, the boiler, the 
heater and the vapor transmission tube being disregarded. 





c d Fig. 4. Pump 2. 
Fig. 3. Pump 1A Fy 2.4 2.5 0.25 cm 
Pump 1B 30 3.5 14 0.15 
Pump 1C 3.5 4.3 0.25 


The pumps which have been tested are listed with their measured data in 
Table I. Some of them are schematically shown in the diagrams (Fig. 3 to 
Fig. 13). Pump 7 (Fig. 9) was constructed in this laboratory by imitating 


imm hole 
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Fig. 5. Pump 3. Fig. 6. Pump 4 (Steel). Fig. 7. Pump 5 (Two stages). 


Gaede and Keesom’s’ high speed pump. Pump 8 (Fig. 10) is Hickman and 
Sanford’s design but with a bigger nozzle chamber and with a removable con- 
denser. Pumps 10, 11, 12 and 13 are new designs of this laboratory.® 


* Hickman and Sanford, J. Phys. Chem. 34, 637 (1930). 
7 Gaede and Keesom, Commun. Leiden 18, No. 195a, (1928). 
8 Ho, R.S.I. 3, 133 (1932). 
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For the low speed pumps, the speeds were measured by a flowmeter, 
which is calibrated for air and other gases, and constantly checked by the 
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_ Fig. 9. Pump 7 (Gaede and 


Fig. 8. Pump 6 (Steel). Keesom’s high speed pump). Fig. 10. Pump 8. 
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Pump 9, 3.6 2.2 0.6cm Fig. 12. Pump 10 (4 nozzles). 
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Fig. 13. Pump 13 (Lauritsen’s 19 nozzle pump). 


constant volume method. For the high speed pumps, the speeds were meas- 
ured by the mercury-pellet method and constantly checked by comparison 
with the known speed of a fore pump. 
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The results in Table I all refer to mercury as a fuel. The speeds were 
specified as (1) at the position along the inlet tube and close to the nozzle 
chamber, (2) constant at pressures from 10~* to 10-5mm Hg, (3) against a fore 
pressure of below 0.1 mm Hg, (4) the water cooling was at room temperature 














TABLE I. 
Pump ‘orm Slit area Watts input Speed Speed 
(for air) factor 
1A Straight nozzle 1.54 cm? 115 1.01 L 0.056 
Small chamber 
1p . 9 1.48 140 >. 0.20 
1c = 2.9 140 7 0.17 
2 Divergent nozzle 1.54 115 1.37 0.086 
Small chamber 
3 - 0.27 110 0.57 0.18 
4 Straight nozzle 30.0 200 6.6 0.019 
Steel 
5 2 stages 1.56 200 0.82 0.045 
Small chamber 
6 Small 1.26 300 1.16 0.079 
Steel 
7 Large 227.0 2200 105 0.040 
Steel 
8 Small : 60 3.3 0.11 
Upward 
9A ° 2.4 60 4.0 0.14 
OB Large 6.4 140 14.5 0.20 
Upward 
10 4 nozzles 9.5 140 28.0 0.26 
Upward 
11 4 nozzles 9.5 140 33.0 0.30 
12 7 nozzles 22.8 200 55 0.21 
13 19 nozzles 155 700 320 0.18 











| 
| 
| 
| 
| 


and the temperature difference between the inlet and the outlet water was 
about 2°C. The slit area of each single nozzle pump was calculated by taking 
the shortest distance from the edge of the nozzle to the condensing wall. The 
watts input includes the very small heat losses by radiation and conduction. 


VI. Discussions AND CONCLUSIONS 
A. Divergent nozzle 


Although the parallel vapor beam cannot possibly be obtained in Craw- 
ford’s pumps? with divergent nozzle as pointed out by Gaede,’ the vapor on 
passing through the divergent nozzle is expanding at the throat and being 
partially rectified by the wall of the nozzle, thus the average upward com- 
ponent of the momentum of the vapor molecules after leaving the nozzle is 
greatly reduced. Consequently, the relative amount of the vapor dispersed 
through the slit into the chamber is diminished. This, combined with the 
effective downward vapor stream, will make the speed factor greater. This is 
why the speeds of pumps 1, and 8 are respectively higher than those of the 
pumps 2 and 9,4, and why the divergent jet is adopted in the design of a high 
speed pump. 


® Crawford, Phys. Rev. 10, 557 (1917). 
10 Gaede, Zeits. f. tech. Physik 4, 348 (1923). 
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B. Chamber resistance 


The chamber resistance which affects mainly the coefficient 8 of the speed 
factor consists of two factors, as mentioned in III, one due to the dimensions 
of the chamber and the other due to the vapor in the chamber. In order to 
minimize this resistance it is necessary to enlarge the space near the slit and 
to keep the chamber cold.. The speed is low if the nozzle is thrust deeply into 
the condensing tube as in the pumps 1,, 2 and 5; while the designs of pumps 
8 and 9 are of low chamber resistance; and their speeds would be still higher if 
they were of “vertically downward form”, because in the “vertically upward 


form” the condensed mercury and the boiler underneath heat the chamber 
considerably. 


C. Slit area 


From the data we can see that the speed factor increases with the increas- 
ing of the width of the slit, if the other parts of the design do not change much. 
This is due partly to the increasing of the coefficient y because the number of 
the gas molecules reflected back through each unit area per unit time is de- 
creased, and partly to the increasing of the coefficient a, since \/a is de- 
creased. Gaede explained this by the term of “extracting power” (Aufnah- 
mefihigkeit) of the vapor stream. Of course, there must be given a maximum 
value to the product ayA by varying A. This is shown partly by pumps 1g, 
1c and 4 which are all of straight nozzle. The slit of pump 4 might be too large 
and that of pump 1g, too small, since their speeds are not as high as they 
should be by comparing their sizes with pump Ic. 


D. Large pumps 


To obtain a high speed we have to enlarge the size of the pump. But the 
high speed cannot be obtained by enlarging all the parts of a small pump in 
proportion. Obviously, the width of the slit cannot be enlarged much if the 
other parts are enlarged in proportion. So the slit area of a large pump is 
usually small in comparison with its size (i.e., the area of the cross-section of 
the condensing tube.) There have been two ways to enlarge the efficient area. 
(i) By using a convergent nozzle with a very narrow throat (Fig. 9) and apply- 
ing a high temperature vapor, Gaede and Keesom were able to get an effective 
vapor stream to operate under a large slit area. (ii) By using several small 
nozzles in place of a single large one and leaving sufficient space between the 
nozzles for the additional path of the gas (Fig. 13), an efficient enlarged slit 
area is obtained. The large speed factor and low power input of the multiple 
nozzle pumps are remarkable. 


VII. REMARKs ON MEASUREMENTS OF SPEED 


Whatever method is used to measure the speed of a diffusion pump the 
pressure of the gas under extraction must be known. In order to get the actual 
speed directly, the pressure should be measured at the inlet mouth of the 
pump. If the pressure gauge is connected at any other point than along the 


inlet mouth a very serious departure from the true speed will sometimes 
occur. 








394 T. L. HO 


(a) Under-value speed 


The speed of a pump can be obtained by measuring the apparent speed 
at the far end of a connecting tube of known dimensions and applying resist- 
ance® equation to calculate the real speed. This method is obviously limited 
by the condition that the connecting tube should be so large that the ideal 
maximum speed, defined by Eq. (1), through an area equal to the area of 
cross section of the tube is much higher than the speed of the pump. Other- 
wise the measurement can never reach the real speed of the pump. The 
author once overlooked the above condition and got a speed of 11 liters in- 
stead of 105 liters per second of the pump 7, using a connecting tube of only 
1.6 cm in diameter. 


(b) Over-value speed 


The departure due to the over-correction is also considerable, if the pres- 
sure gauge is directly connected to the nozzle chamber. Gaede stated, for his 
three-stage steel pump, that at a point adjacent to the slit a speed of 60 liters 
per second was obtained, but at the wide mouth of the inlet tube the added 
resistance (i.e., chamber resistance) reduced the result to 15 liters per second. 
From his statement it is seen that 15 liters per second is the real speed of his 
pump; and 60 liters per second is the over-value speed, but it would have been 
the true speed if the nozzle chamber and the inlet tube had been enlarged and 
cooled. 


(c) **Pseudospeed”’ 


As in above, if the pressure gauge is directly connected to the nozzle 
chamber of a pump, the inlet of which is not along the same axis as the 
chamber, the pressure thus measured is not on the main path of the gas and 
will render the value of the speed meaningless. For, not only can the pressure 
not be the same throughout the nozzle chamber because of the chamber re- 
sistance, but also the pressure distribution cannot be symmetrical with re- 
spect to the axis of the chamber (or of the nozzle), especially for the pump of 

















TABLE II. 
Speed (for air) Pseudo-speed (for air) 
Water cooling Still air Water cooling Still air 
Pump : ——___—_—_—_————-cooling (cop- |——-—-—— ———~— cooling (cop- 
With cooling Without cool- per wire With cooling Without cool- per wire 
around ing around projection) around ing around projection) 
chamber chamber chamber chamber 
& a.514, 2.8 0.17 3.8L 3.2 0.50 
OR 14.5 10.0 25.0 28.0 
10 28.0 15.0 64 70 











“vertically upward form” with a small and hot chamber. The pseudospeeds 
of pumps 8, 9g and 10, measured by connecting the pressure gauge at a point 
opposite to the inlet tube, are shown in Table II. It is interesting to note that 
the pseudospeeds of pumps 9, and 10 are several times larger than their cor- 
responding true speeds and decreases with the decreasing of the temperature 
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of the chamber. This is because the hot vapor raised from the lower part of 
the chamber prevents the gas passing from the inlet tube to the opposite 
side (where the gauge is connected) so that the partial pressure of the gas 
(which is measured by the gauge) on that side is low. It can of course be in- 
creased by cooling the chamber. The pseudospeed of pump 8, because of the 
large chamber and its great distance (8 cm) from the boiler, is nearly the true 
speed when the pump is properly cooled. 

Recently, Hickman and Sanford’ measured the speeds of their pumps by 
connecting the McLeod gauge to the chamber at a point far from the inlet 
mouth of the pump (seen from their schematical diagram of the apparatus 
and their pictures of the pumps). Some of their speed-pressure (fore) curves 
(for instance, graph 1 and graph 6 in their paper) show that the speed of a 
pump by “still air cooling” is equal to or higher than that by “fan cooling”, 
the other conditions being the same. This is inexplicable and in contradiction 
to Crawford’s experimental results,® and also to one of the general experi- 
mental results (i.e., (e) of section IV). Further, as the author tested, for the 
mercury filled pump the speed is very low if it is only cooled by the still air 
(with copper wire projections). The speeds Hickman and Sanford measured 
are probably pseudospeeds which are higher than the corresponding true 
speeds. 

The pseudospeeds shown in Table II are, as true speeds, independent 
of the pressure of the extracted gas. This shows that the pressure distribution 
inside the nozzle chamber is the same for different pressures. 

In conclusion, the author wishes to express his thanks to Professor A. 
Goetz and Professor W. R. Smythe for their valuable criticisms, and to 
Professor C. C. Lauritsen for his helpful suggestions and for the use of his 
pumps. 
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An A.C. Wheatstone Bridge for Audio and 
Radio-Frequency Measurements 


By A. Hemincway and J. F. McCLENDoN 


Laboratory of Physiological Chemistry, University of Minnesota, Minneapolis 
(Received March 12, 1932) 


N MANY biological problems it is necessary to measure the conductivity 

and capacity of suspensions, solutions, and tissues with varying fre- 
quencies of the alternating current. The method of electrical conductivity 
and capacity is perhaps the best method of investigating the phenomenon 
of polarization, which, according to Nernst! is the fundamental principle of 
electrical excitation of tissues. This method has been used by Gildemeister,? 
Hozawa,’ Lullies,t and others. According to Osterhout® and Cowdry,® the 
electrical conductivity method is the best quantitative method of following 
permeability changes of tissues. This is the method employed by McClen- 
don,’ Brooks,* Osterhout.’ Fricke,!° and McClendon!!” have used the 
method of electrical conductivity measurements to study the physical struc- 
ture of blood cells and Fricke and Stewart have used this method to meas- 
ure the volume of the disperse phase in biological suspensions. In a quantita- 
tive investigation of heat dosage in diathermy Hemingway measured the 
radio-frequency resistance of the human body" and living tissues."® 

These few examples will illustrate the various biological uses of a Wheat- 
stone bridge, which will be suitable for the measurement of resistance and 
capacity at various frequencies. The building of a suitable bridge presents 
considerable difficulty if measurements are to be made at frequencies higher 
than those of the audio-range. Philippson,'7'* and Cole'® have measured by 


1W. Nernst, Pfluger’s Arch. 122, 275 (1908). 

2M. Gildemeister, Pfluger’s Arch. 219, 89 (1928). 

%S. Hozawa, Jour. of Biophysics 1, 185 (1925). 

* Lullies, Pfluger’s Arch. 221, 298 (1928). 

5 W. J. V. Osterhout, Jour. Biol. Chem. 36, 485 (1918). 

6 E. V. Cowdry, Cytology, University of Chicago Press (1924). 

7 J. F. McClendon, Am. Jour. Physiol. 29, 302 (1912). 

8S. C. Brooks, Jour. Gen. Physiol. 7, 349 (1925). 

* W. J. V. Osterhout, Injury, Recovery and Death in Relation to Conductivity and Per- 
meability, Philadelphia (1922). 

10 H. Fricke, Phys. Rev. 26, 682 (1925). 

J. F. McClendon, Protoplasma 3, 7 (1927). 

2 J. F. McClendon, Protoplasma 3, 71 (1927). 

‘8H. Fricke, Jour. Gen. Physiol. 6, 741 (1924). 

4G, N. Stewart, Am. Jour. Physiol. 90, 194 (1929). 

1 A, Hemingway, Radiology 14, 99 (1930). 

16 A. Hemingway and-D. C. Collins, Am. Jour. Physiol. 94, 338 (1932). 

17 M. Philippson, Compt. rend. Soc. Biol. 83, 1399 (1920). 

‘8 M. Philippson, Arch. int. Physiol. 18, 161 (1921). 

1? K. S. Cole, Jour. Gen. Physiol. 12, 29 (1928). 
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resonance methods the impedance at various frequencies. The bridge method, 
however, enables the experimenter to separate the resistance factor of the im- 
pedance from the reactance factor, which is of importance in biological work. 
In the past few years” we have been building and perfecting a Wheatstone 
bridge suitable for the type of measurement mentioned and it is this bridge 
which is described below. 


THE BRIDGE CIRCUIT 


The circuit used is similar to that devised by Stratton,?° and utilized by 
Remington.” [t is shown in Fig. 1. The four arms of the bridge contain: 1 and 
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Circuit Diagram of Wheatstone Bridge 
Fig. 1. 


2, equal ratio arms; X, the unknown whose resistance and capacity is de- 
sired; and C, R, r, a balancing arrangement as shown of two variable decade 
resistors and a decade capacity. The bridge is also designed to allow an in- 
ductance to be placed in series with the unknown in arm X, an arrangement 
used by many investigators. 


20 J. A. Stratton, Jour. Opt. Soc. Am. 13, 471 (1926). 
21 R, E. Remington, Protoplasma 5, 339 (1928). 
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For measurements in the audio-frequency range a Vreeland oscillator is 
used as a source of alternating current. The characteristics of the Vreeland 
have been investigated by Taylor and Acree,” Morgan and Lammert”®’ and 
others and found to be satisfactory for a. c. measurements due to the minimum 
of alternating current harmonics. When the Vreeland oscillator is connected to 
the bridge, an audio-frequency amplifier with ear phones is used in the bridge 
circuit for balancing. 

For radio-frequency measurements a heterodyne method is used. The 
main oscillator and the heterodyne oscillator produce alternating currents 
whose frequencies differ by approximately 1000 cycles per second, and these 
two frequencies are impressed on the detector. The plate current of the de- 
tector passes through the primary of the first audio-frequency transformer of 
the audio-frequency amplifier, which amplifies the beat note of 1000 cycles 
per second. 

Hitherto for our work it has been sufficient to make measurements of re- 
sistance and capacity throughout the audio-range of frequencies, and to make 
a radio-frequency measurement at a single radio-frequency of one million 
cycles per second. It is our intention, however, in the future to install the 
proper inductances and capacities in the oscillators and detector circuit for 
measurements at intermediate frequencies. 


THE OSCILLATORS, DETECTOR AND AMPLIFIER 


For oscillators the two tube Eccles-Jordan type recommended by Strat- 
ton®® has been found to be very satisfactory. The main oscillator has two 
UX112A tubes with 150 volts on the plates while the heterodyne oscillator 
has 50 volts on the plates of two UV201A tubes. The detector has a UV201A 
tube with a plate voltage of 25 volts. An important and necessary feature of 
these circuits is the inductances which are wound in a toroidal form to reduce 
the external field. The condenser of the main oscillating circuit is fixed, while 
those in the detector and heterodyne oscillator are variable for tuning. In the 
oscillating circuit of the main oscillator is placed a Weston thermocouple 
milliameter which indicates whether oscillation is taking place and whose cur- 
rent value indicates whether or not the batteries need charging. The amplifier 
is of the usual type with two stages of amplification. UV201A tubes are used 
with 50 volts on the plates and transformer coupling. 


RESISTORS, CONDENSERS, AND RATIO ARMS 


Two decade resistors are used for balancing, one for radio-frequency and 
the other for audio-frequency, and are connected as shown in Fig. 1. For radio- 
frequency measurements the 1, 10, and 100 ohm decades are wire woven 
coils, with a modified Ayrton-Perry winding and constructed especially for 
our purpose by members of the Radio Engineering Department of the Uni- 
verstiy of Minnesota. The 1000 and 10,000 ohm units are “El Menco” carbon 
coated cylindrical resistors. These are calibrated with 1000 cycle alternating 


2 W. A. Taylor and S. F. Acree, Jour. Am. Chem. Soc. 38, 2396 (1916). 
* L. R. Morgan and O. M. Lammert, Jour. Am. Chem. Soc. 48, 1220 (1926). 
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current when necessary and are kept at fairly constant temperature. For the 
lower frequency measurements of the audio-range a General Radio type 102 
decade resistor is used. The condensers used are Leeds and Northrup preci- 
sion decade condensers with a variable air condenser for capacities less than 
0.001 microfarad. 

Ratio arms consist of two “El Menco” units of 100, 1000, or 10,000 ohms 
for high frequency or two wire resistors wound on long cards with the modi- 
fied Ayrton-Perry method for audio-frequencies. The two coils are connected 
in series with three contacts, one from each end of the pair and one from the 
center, which project downwards into cups of mercury. This method allows of 
easy interchange of the arms and convenient reversal. The use of equal fixed 
ratio arms eliminates the difficulties due to drum wound slide wires, men- 
tioned by Jones and Joseph,** Morgan and Lammert,” Schlesinger and Reed,** 
and Parker.” 

SHIELDING 


Proper shielding of the bridge for accurate measurements is essent 
especially with radio-frequency currents. This matter has occupied the 
tention of engineers and physicists in the past few years and several method 


are now in general use. One is the older method of Campbell’? recently adop- 
ted by Shackelton,?* Ferguson” and Shedlovsky,*’ and is particularly useful 
in transmission line work. In this method the individual bridge units are in- 
closed in a “nest” fashion and the ground currents brought to one or more 
points of the bridge. This method requires a considerable number of metallic 
shielding devices which are unnecessary for the type of bridge described in 
this paper. In the bridge which we have designed the ground leakage has been 
reduced to a minimum by using spacious metallic containers in which the 
parts of the bridge are separated from the sheet metal ground containers with 
a wide air space. This decreases the capacitative short circuiting of a resistor 
by a metallic conductor in close proximity to it. Galvanized sheet iron is used 
for shielding. This was chosen, rather than copper which was used tempor- 
arily, since measurements were to be taken over a wide range of frequencies 
and the slight disadvantage due to the lower resistivity of the sheet iron 
would be compensated by the ferromagnetic properties, which would be 
effective at the lower frequencies.* 


THE WAGNER GROUND 


Since Wagner” first proposed his now famous earthing device, by which 
the leak to ground is made equal on both sides of the bridge, several types 


4 G, Jones and R. C. Joseph, Jour. Am. Chem. Soc. 50, 1049 (1929). 

*% H. I. Schlesinger and F. H. Reed, Jour. Am. Chem. Soc. 41, 1727 (1919). 
% H.C. Parker, Jour. Am, Chem. Soc. 45, 1366 (1923). 

27 G. A. Campbell, Elec. World 43, 647 (1904). 

28 W. J. Shackelton, Trans. A.I.E.E. 45, 1266 (1927). 

29 J. E. Ferguson, Trans. A.I.E.E. 48, 517 (1929). 

30 T, Shedlovsky, Jour. Am. Chem. Soc. 52, 1793 (1930). 

st J. H. Morecroft and A. Turner, Jour. Inst. Rad. Eng. 13, 477 (1925). 

#2 K. W. Wagner, Elektrotech. Zeits. 32, 1001 (1911). 
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have been used, among which are those of Jones and Joseph,” Schlesinger and 
Reed*® and Shedlovsky.*® The one described here is that used by Stratton” 
and consists of two air condensers connected to each of the input leads and to 
ground, with a variable decade resistance of from 1000 to 110,000 ohms, 
which can be switched to either input terminal by a d. p. d. t. switch. This 
arrangement can be adjusted very easily and has proved to be very satis- 
factory. Only by a careful setting of the Wagner ground can consistent re- 
sults be obtained with radio-frequency alternating current. 


THE BRIDGE ASSEMBLY 


The exterior of the assembled bridge is shown in Fig. 2. There are two 
large sections joined by a cylindrical tube. Within the smaller section is the 
main high-frequency oscillator and its batteries. Within the cylindrical tube 
joining the two sections of the bridge are the two leads from the oscillator to 
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Exterior of the Wheatstone Bridge 








Fig. 2. 


the bridge proper. The Vreeland oscillator which furnishes the audio-fre- 
quency current is not shown. 

The larger section shown in Fig. 2 has two stories. The upper contains the 
four arms of the bridge, that is the resistors, capacities, the Wagner ground 
and the constant temperature air bath with the thermostat. The lower story 
is divided into six compartments, three of which are on the side facing the 
main oscillator and these three compartments contain batteries, each a 6 volt 
A battery and B batteries. On the opposite side from the main oscillator of 
this larger section and in the lower story are three compartments containing 
the heterodyne oscillator, the detector and the audio-frequency amplifier. 
Plans showing the arrangements of these two stories are shown in Figs. 3 and 
4 and 5. All the apparatus of the upper story is supported on glass rods whose 
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ends rest on small metallic brackets connected to the sheet iron shield. Be- 
tween the compartments the wires pass through holes two inches in diameter. 
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Diagram of Upper Story of Wheatstone Bridge 
Fig. 3. 


Ebonite or glass rods project from the sliding contacts of the decade dials up- 
wards through small holes in the top of the bridge. 
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oscillator! 


Diagram of Upper Story of Wheatstone Bridge 
Fig. 4. 


Figs. 3 and 4 show the arrangement of the upper story. Fig. 3 is the sche- 
matic wiring diagram and Fig. 4 the actual arrangement. Fig. 4 is shown as 
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viewed from above with the roof removed. A is the low frequency connection; 
the switch in this position leads to the amplifier directly. In the position D the 
connection is to the detector. Both the amplifier and detector are in the lower 
story as in Fig. 4. A detailed diagram of the arrangement of the main oscil- 
lator, which is the smaller of the two sections, is not given. 


ACCURACY 
At radio-frequency the accuracy is about 1 percent. To test the bridge 
manganin wire standards are used. These consist of a single strand of thin 
manganin wire which runs back and forth between pegs set on opposite edges 


of a large hardwood board about 1 inch apart. Wound in this way the dis- 
tributed capacity and inductance are a minimum and the skin effect is small 
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Fig. 5. 


due to the relatively high resistivity of the wire and its small diameter. With 
these standards the value of the radio-frequency resistance equals that of the 
1000 cycle or direct current resistance to within 1 percent. The audio-fre- 
quency measurements are accurate to 0.1 percent. 

We wish to take this occasion to thank Professor A. Zeleny of the Depart- 
ment of Physics and Professor H. E. Hartig of the Department of Electrical 
Engineering of the University of Minnesota for the use of apparatus for col- 
laboration purposes and for their asistance and cooperation in this work. 
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Free Rotation in Molecules 


HE rotation of molecular groups (radi- 

cals) within complex molecules, about 
the bonds which hold them together is a phe- 
nomenon to which considerable interest is at- 
tached. In forthcoming issue of the Physical 
Review a paper! by H. Nielson gives an in- 
teresting discussion of the quantum theory of 
an effect of this kind. In a previous paper? in 
the same journal, L. Pauling has considered 
the theory for an analogous case involving the 
vibration and rotation of molecules in a 
crystal lattice. Nielson takes as a model a 
body dynamically symmetrical about an axis 
of figure but divided into two parts by a plane 
perpendicular to the figure axis. The two parts 
thus formed are assumed to be capable of ro- 
tating with respect to each other about the 
figure axis but subject to a restoring torque 
tending to keep them in an equilibrium posi- 
tion. A rough mechanical model would be two 
identical right cylinders placed end to end and 
joined by a wire passed through them along 
their common axis, supposing that the wire 
could twist but not bend, the wire playing the 
roéle of the chemical bond. There is one im- 
portant difference between this mechanical 
model and the actual molecule, namely that in 
the former there is but one definite equilib- 
rium position being that in which the wire is 
not twisted at all, for which the restoring 
torque vanishes, but in the latter there may 
be two or more equilibrium positions depend- 
ing upon the symmetry of the arrangement of 


the atoms within the molecule. In order to 
take this into account Nielson supposes that 
the restoring torque can be obtained from a 
potential of the form V(¢’)~(1—cos md¢’) 
where ¢’ is the angle of relative rotation of the 
two parts and m is an integer. In the range 
0 <¢’ S2r this function vanishes for 2m values 
of ¢’, there being m such values in the range 
0<¢’<z and m further values symmetrically 
placed in the range t<¢’ <2r. 

When Schrédinger’s equation is written 
down for such a system the part depending on 
the angle #’ can easily be separated out and 
is found to be the standard equation for the 
Mathieu functions familiar in mathematical 
analysis. If one restricts @’ to have values 
only in the neighborhood of one of the equi- 
librium positions then the equation reduces to 
that for the harmonic oscillator, in agreement 
with the fact that for small oscillations the 
vibrations of a torsion pendulum are simple 
harmonic. The complete solution of the equa- 
tion presents difficulties arising from the 
periodic character of the potential function 
which is reflected in the periodic character of 
the solutions of the Mathieu equation. The 
positions of the energy levels and the transi- 
tion probabilities are also discussed by the 
author. 


E. HILi 


1H. Nielson, Phys. Rev. 40, 445 (1932). 


2 L. Pauling, Phys. Rev. 36, 460 (1932). 


The Neutron Again 


N VIEW of the importance and interest in 

the subject it seems suitable to amplify a 
little the few remarks made in these columns 
last month on the subject of the neutron.! In 
particular we must mention the important ex- 
periments of Mme. I. Curie and M. F. Joliot 
which form one of the strongest links in the 
present chain of evidence pointing towards 


1 Editor’s Column, Physics 2, 310 (1932). 


the neutron.?:* Since the appearance of Chad- 
wick’s note‘ they have presented further evi- 


21. Curie, Comptes Rendus 193, 1412 
(1931); F. Joliot, Comptes Rendus 193, 1415 
(1931). 

37. Curie and F. Joliot, Comptes Rendus 
194, 273 and 708 (1932). 

‘J. Chadwick, Nature Feb. 27 (1932) p. 
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dence of a very convincing nature, convincing 
at least of the peculiarity of this new radia- 
tion.» They have made a very careful study of 
the interesting way in which this radiation is 
absorbed by substances containing protons or 
other light nuclei such as helium. The radia- 
tion to which we refer is that produced from 
beryllium atoms when they are bombarded by 
the a-rays from polonium. The anomalous 
manner in which it is absorbed on passing 
through paraffin, cellophane, etc. even though 
its absorption coefficient in lead is very small 
makes it quite clear that this radiation is not 
of an electromagnetic nature. Both theory and 
experiment agree that the ability of a charged 
particle to take up energy from electromag- 
netic radiation varies in a strongly inverse 
manner with the mass of the particle. Thus in 
the absorption of x-rays by matter, the effect 
is produced almost entirely by the electrons 
within the atoms, the nuclei themselves pro- 
ducing only a small effect. But in the present 
case the reverse is true, and one proton is 
found to be as efficient as four electrons in 
absorbing the radiation. 

In a recent note evidence of a visual nature 
has been presented which is quite appealing in 
its directness.’ A large Wilson cloud-track 
chamber is placed in the direct path of the 
radiation, and within it is placed a mica screen 
covered with paraffin on which the radiation 
falls after being filtered through 2 cm of lead. 
When the chamber is operated tracks are ob- 
served due to protons knocked out of the 
screen by the radiation, these tracks emanat- 
ing in a bundle from the screen. But in addi- 
tion to these a number of other tracks have 
been observed which do not start from the 
screen but from some point within the body 
of the chamber itself, that is, from a point 
within the “water vapor” region of the cham- 
ber. These tracks appear to be due to nuclei 
or protons which have experienced a collision 
with another particle which itself does an in- 
appreciable amount of ionizing and so leaves 
no track in the cloud chamber. Photographs 
of such tracks have been made and cuts of 
them are shown in the publication. This all 
agrees very well, at least qualitatively, with 
the interpretation of the radiation as being 
composed of neutrons. Being uncharged a 
neutron would have a very restricted field of 
force, and a rather close approach would be 
required for a collision with an atom to be 
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effective in producing ionization. The neutron 
might travel quite a distance bet ween ionizing 
collisions. 

Similar experiments have been successfully 
carried out by P. Auger*® who has studied the 
paths followed by the ejected protons in a 
magnetic field. He estimates that the protons 
have energies from 5X10* up to several 
million electron-volts. 

In the previous note! we neglected to men- 
tion one important bit of evidence adduced by 
Chadwick‘ and Webster;’ i.e., the experimen- 
tal observation that the rays emitted from the 
bombarded beryllium atoms in the same di- 
rection as the incident a-particles travel is 
harder (more penetrating) than the radiation 
emitted in the reverse direction. This is easily 
understandable on the hypothesis that the 
radiation is composed of neutrons behaving 
like particles of unit atomic weight, as the 
conservation laws for energy and linear mo- 
mentum show that the neutrons projected in 
the forward direction by the a-particles would 
have about 50 percent more energy than those 
ejected the backward direction. On the other 
hand, one might think that even if the radia- 
tion were made up of photons a similar result 
might hold since Chadwick’s computations 
show that the photon would have an energy 
of around 14 10° electron-volts, and such a 
photon would have a mass equal to about 28 
times the rest mass of an electron. As it moves 
with the velocity of light it carries away an 
appreciable fraction of the momentum of the 
incident a-particle. If the momentum of the 
photon varied with its direction of projection, 
then its frequency would do likewise, and the 
radiation emitted in different directions with 
respect to the direction of the incident a-par- 
ticles would have different absorption coef- 
ficients. However a direct computation of the 
difference in frequency for the forward and 
backward directions of emission of the photon 
shows that it is negligible under the conditions 
of the experiments. The evidence thus favors 
the neutron hypothesis. 

E. HILu 


57. Curie and F. Joliot, Comptes Rendus 
194, 876 (1932). 

6 P. Auger, Comptes Rendus 194, 877 
(1932). 

7H. C. Webster, Nature, March 12 (1932) 
p. 402. 
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Cosmic Rays: What Are They? 


CHARGED electroscope placed in air 

will not keep its charge no matter how 
carefully the instrument is insulated. The 
rate of “leak” is small but too large to be 
simply due to the movement of the charge 
along the surface of the insulator. The fact 
that the rate of loss of charge does not depend 
on the dimensions of the insulator and other 
tests prove conclusively that the air surround- 
ing the charged body is ionized. These air ions 
then discharge the electroscope. A part of the 
ionization in air is produced by the relatively 
small amount of radioactive material in the 
air, the underlying soil and rocks and in the 
material of the electroscope itself. But after 
all known sources of ionization are taken into 
account a residual ionization remains. “Cos- 
mic rays are by definition, the cause of that 
residue of the ionization.”! 

Physicists in various parts of the world 
have been carrying on experiments to deter- 
mine the nature of the ionizing agent that 
produces this residual ionization. It was very 
soon discovered that the ionizing agent what- 
ever its nature, had an extremely high pene- 
trating power. Regener? reports detectable 
ionization 186 meters and 230 meters below 
water. The most penetrating of the known 
gamma-rays spontaneously emitted from ra- 
dioactive materials could hardly be detected 
through 20 cm of lead. If one assumes that the 
ionizing agent observed by Regener is a short 
electromagnetic wave (or corpuscle, as you 
please) whose absorption depends on the mass 
through which it travels, 200 cm of lead would 
be required to reduce the intensity to the 
minimum value which he observed. 

To account for this penetrating power, 
Millikan and his school have in the past as- 
sumed that the cosmic rays were radiations 
similar in nature to gamma-rays except much 
harder, that is, of much higher frequency. 
Evidence supporting this hypothesis is the 
constancy of the radiation at all parts of the 
earth’s surface. Millikan and his co-workers 
have made measurements in California, in the 
Bolivian High Andes near the equator and at 
Churchill, Manitoba, from the 
North magnetic pole. If the rays were elec- 
trons their distribution at the earth’s surface 
would be modified by the earth’s magnetic 
field. To test further the constancy of the in- 
coming radiation over various parts of the 
earth, Dr. A. H. Compton and associates are 


730° miles 


now taking readings at widely separated 
points, both at high and low altitudes and at 
high and low latitudes in the northern and 
southern hemispheres. 

We have tacitly assumed that the rays 
come from outer space into the earth. Obser- 
vations with an electroscope on the surface 
of the earth, however, do not determine the 
direction of travel of the rays. If the source of 
the rays is on the earth’s surface one would 
expect the rate at which ions are produced in 
an electroscope to be decreased as one as- 
cended upwards, leaving an absorbing blanket 
of air between the instrument and the earth. 
When balloons carrying instruments are sent 
up 6 to 10 miles it is found that the ionization 
increases. At sea levels between one and two 
ions are on the average produced per cc per 
sec. in air under standard Pic- 
card, in his famous balloon trip found about 
200 per cc per sec. in air (at N.P.T.) at an 
altitude of 10 miles. Such evidence with that 
taken below the surface of lakes, rules out the 


conditions. 


possibility of a terrestrial source for the rays. 

Millikan interprets the manner in which the 
ionization varies with height to be evidence in 
support of the photon nature of the radiation. 
The ionization produced when photons pass 
through a gas is mainly produced by the sec- 
ondary electrons which are ejected from the 
gas molecules. At very high frequencies of the 
photons, theory and experiments indicate that 
a large fraction of the electrons ejected from 
the gas atoms will be “Compton-effect” elec- 
trons and many of them will follow very 
closely the original direction of the photon 
beam. These electrons will have a definite 
range in a gas, and until one descends below 
the “top of the atmosphere” to a distance at 
which the number of the secondary electrons 
are removed from the path or its immediate 
neighborhood is equal to the number of them 
that are produced, the ionization per unit path 
length will increase. Below such a point the 
total ionization produced would be decreased 
since the number of primary photons would 
be reduced. At the altitude at which the 

! Darrow, Bell System Tech. Jour. Jan. 
1932, gives a complete review and discussion 
of this subject, covering the literature to a 
late date in 1931. 

2 Regener, Die Naturwissenschaften 15, 


183 (1929). Zeits. f. Physik 74, 433 (1932). 
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beam of photons produces a maximum ioniza- 
tion in a closed vessel, the number of ionizing 
sources removed from the beam is equal to the 
number produced in it. At that point, the 
beam is said to be in equilibrium with its sec- 
ondaries. A variation with height in a.manner 
consistent with the photon hypothesis has 
been found in the balloon experiments. The 
number of ions produced increased to a maxi- 
mum at a given height and then decreased at 
higher altitudes. 

With the assurance that the cosmic rays 
were electromagnetic waves, Millikan has 
used various absorption coefficients to get an 
explanation of his curves giving the variation 
of ionization with depth below water. To get 
a calculated curve to fit the empirical curve, 
he assumes three or more groups of waves 
having different absorption coefficients. By 
means of formulae relating the absorption 
coefficient to frequency, the frequencies of 
various cosmic rays are computed and by the 
use of E =hv, their energies are given. To find 
a source for the energy of the magnitude 
(2108 electron volts) given by these calcu- 
lations presented some difficulty. Millikan 
suggested that “in the depths of space” the 
celestially common elements, helium, oxygen 
and silicon were created by the union of the 
proper number of protons and electrons. In 
this process some mass is lost as is indicated 
by Aston’s curves of atomic masses. If one 
assumes that the disappearance of a mass M 
is accompanied by the appearance of radiation 
whose energy is given by MC?=E where 
E=hy is the energy and C=3 X10! cm/sec., a 
correlation between the energies of the cosmic 
rays and the suggested atom-building proces- 
ses is obtained. This suggestion has led to 
numerous speculative discussions concerning 
the “cycle of the Universe”. It allows one to 
assume the mass of the stars to be continually 
changing into radiation without necessarily 
assuming that the final state of the Universe 
is one of radiation existing alone. Since the 
cosmic rays do not vary with the position of 
the observer with respect to the Milky Way 
or any other celestial body, it is assumed that 
the atom-building processes take place in free 
space away from the stars. And the suggestion 
is made that possibly a process the reverse of 
that taking place in the stars, is going on in 
the depths of space, that is radiant energy 
condensing into atoms.* 

Regener? concludes that the most penetrat- 
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ing of the radiation he detected arises from the 
annihilation of an entire He atom and the 
creation of electromagnetic radiation. The 
wave-length of such radiation is about three 
millionths of an angstrom unit. 

The above conclusions are all based on the 
results obtained in closed vessels filled with 
gases at atmospheric or higher pressures. The 
total ionization in the vessel is determined and 
no information given concerning the details of 
the production of this ionization. By use of the 
Wilson cloud expansion apparatus it is possi- 
ble to detect visually the path of individual 
ionizing agents, in particular alpha particles, 
swift electrons and protons. Skobelzyn several 
years ago observed paths in an expansion 
chamber in which gamma-rays were liberating 
secondary electrons. A magnetic field was ap- 
plied to bend the electrons in spirals and from 
the radius of curvature of the spiral, the 
speed of the electrons could be measured. A 
magnetic field was applied to bend the elec- 
trons in spirals and from the radius of cur- 
vature of the spiral, the speed of the electrons 
could be measured. The gamma rays were 
those from known radioactive materials. Some 
of the paths were far too straight to be ac- 
counted for by electrons produced by the gam- 
ma rays. This observation led to similar ex- 
periments to study the cosmic rays, since it 
was assumed that these straight paths were 
caused by the cosmic rays themselves or by 
secondary rays which in turn were produced 
by the cosmic rays. 

Mott-Smith and Locher* combined a Wil- 
son cloud expansion apparatus with a Geiger- 
Miiller counter. Bothe and Kolhérster® had 
shown that a Geiger-Miiller counter would 
respond to individual ionizing particles (what- 
ever their nature). These physicists had placed 
one counter above another and obtained 
numerous simultaneous discharges in them. 
Mott-Smith and Locher desired to learn 
whether or not the agents causing the simul- 
taneous discharges in the counters were the 
same as those producing the visible paths in 
the apparatus of Bothe and Kolhérster. They 
placed the Wilson expansion chamber below 


3 Millikan and Cameron, Phys. Rev. 32, 
533 (1928). 

* Mott-Smith and Locher, Phys. Rev. 38, 
1399 (1931). 

5 Bothe and Kolhdérster, Zeit. f. Physik 56, 
751 (1929). 
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the Geiger-Miiller counters and concluded 
that visible tracks were produced by the same 
agents that produced the simultaneous dis- 
charges in the Geiger-M iiller counters. Such 
visible straight line tracks have previously 
been detected only when produced by ma- 
terial particles; protons, alpha-particles or 
electrons. 

Does this then mean that cosmic rays are 
charged material particles? If so magnetic 
fields of very high strength may be used to 
measure their velocity and energy. Curtiss 
earlier had obtained results which indicated 
that the number of coincidences of discharges 
in two Geiger-Miiller counters was reduced 
when a transverse magnetic field was placed 
between them. Mott-Smith® used three 
Geiger-M iiller counters placed in a vertical 
line and observed triple coincidences. The 
first two counters served as a “slit system” 
defining a beam which passed through the 
third. Between the second and third counters 
a block of iron is placed which can be mag- 
netized. Observations were then made to de- 
termine the change in the number of triple 
coincidences when the iron was magnetized. 
The results were essentially negative, no 
change in the number of coincidences being 
observed. 

This result may be interpreted as indicating 
that the particles have an energy of 210° 
electron-volts if electrons or half that amount 
if protons energies much higher than involved 
in the calculations on the assumption that the 
rays were electromagnetic waves. Work by 
Rossi indicates even higher energies. One way 
of accounting for the straight paths is to as- 
sume that the particles are neutron particles 
having zero charge and probably the same 
mass as the proton. Having zero charge they 
are not deflected by the magnetic field. In all 
the work with iron it is assumed that the field 
effective in deflecting a charged particle is 
equal to the magnetic induction B. If the 
effective field should be the magnetizing field 
H no deflection would be expected. A test of 
whether B or H is to be used in such experi- 
ments might be made by sending swift beta- 
particles through a thin sheet of magnetized 
iron. 

The visible paths and the discharges of the 
counters may be caused by particles of sec- 
ondary origin, produced by a primary ultra- 


& Mott-Smith, Phys. Rev. 39, 391 (1932). 
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gamma ray. Such an explanation has some 
support in that some of the paths in the cloud 
expansion apparatus seem to be produced 
simultaneously and to have a common start- 
ing point. Even if these particles were neu- 
trons they may have been produced by high 
energy photons. If neutrons are produced by 
bombarding atoms with alpha-particles (see 
preceding note in these columns) may not 
ultra-gamma-rays be able to do the same? 


To accept the secondary origin for these 
rays one must reconcile the fact that the 
curves showing the reduction of number of the 
visible paths and coincident discharges of the 
G—M counters have the same form as those 
giving the ionization in closed vessels versus 
depth in water. Unless the range of the second- 
ary rays produced by the ultra-gamma-rays is 
very short a different form of curve would be 
expected. The range of the particles however 
is not short. If the observed rays are second- 
ary in nature one would expect their average 
energy to be less than that of the initial pho- 
ton. But the energies mentioned just above 
are much higher than those indicated by the 
gamma ray absorption calculations. 

In a short article published in the May first 
issue of the Physical Review Millikan and 
Anderson’ describe experiments in which 
tracks 15 cm long are produced by means of a 
magnetic field of 17,000 to 20,000 gauss in air. 
Out of 1000 exposures, 34 show measurable 
cosmic-ray tracks and on 24 of these photo- 
graphs tracks of positive particles are found. 
On 6 of the photographs associated tracks are 
found, four of these plates showing tracks of 
one positive and one negative particle. The 
energies calculated on the basis of the particles 
being electron and proton are found to range 
from 16 to 500 million electron-volts with 
possibly a few of them having energy as high 
as 1000 million electron-volts. These authors 
conclude that the primary cosmic rays are not 
neutrons but are photons with energies corre- 
sponding to the atom-building of oxygen, sili- 
con and iron. They believe that most of the 
tracks observed are produced by the bombard- 
ment of the nuclei of the air molecules by 
these photons. 


The data thus are contradictory, at least 
the interpretations are, and the question as 


7 Millikan and Anderson, Phys. Rev. 40, 
325 (1932). 





408 


to the nature, energy and source of that 
which has been called cosmic rays becomes 
more puzzling than 


ever. Perhaps these 


rays are combinations of photons electrons 
neutrons and possibly protons and the vari- 
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